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ABSTRACT 

A thorough understanding and a well developed fabrication 
procedure of MIS structures are the prerequisite for CCD appli- 
cations. The object of this thesis is to study the narrow gap 
semiconductor MIS and investigate its feasibility for IR-CCI 
applications. Two MIS studies were made. 

First, MIS of five lead-tin semiconductors were fabricated 

O 

using E-gun evaporated 100-450A thick AI 2 O 3 or Si 02 layers as 
insulators. C-V measurements Indicated that these MIS behave 
qualitatively like that of Si-MOS. Accumulation, depletion 
and inversion layers were controlled by the gate voltage. How- 
ever, comparisons of measured C-V with theoretical calcula- 
tions did not yield quantitative agreement. 

Second, MIS of 0.05^-cm p-type and 40J^-cra n-type Ge were 
also studied. C-V and C-t measurements indicated standard MOS 
behavior although some small hysteresis was found. Analysis 
based on C-V data showed that the flatband voltage was appro- 
ximately - 1.2 volts' and interface state density was on the order 
of 10 /cm^ . 

In addition, effects of electron bombardment simulating 
the space environment around Jupiter on a n-channel depletion 
MOSFET were studied. The negative threshold voltage was de- 
creased and at a total dose of 2X10 e/cm^ , it became posi- 
tive making the MOSFET an enhancement type of very poor 
quality . 
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I. INTRODUCTION 



The charge-coupled device (CCD) concept was first intro- 
duced by Boyle and Smith of Bell Telephone Laboratory 
around March, 1970 [2]. A new device principle was proper 
in which the signal was represented by minority carriers 
in the depletion layer of a metal-oxide semiconductor 
structure. It is a new device concept but uses Si-MOS 
technology already in existence. 

The development in the past three years was indeed 
amazing and CCD has demonstrated its major impact on three- 
electronic areas: imaging, analog delay and memory. The 
influence is different in these areas, however. 

(1) For imaging, the impact is major mainly because 
the competing imaging technologies all used, in one way or 
another, the relatively more complicated electron-beam 
scanning mechanism and high-vacuum technology of the con- 
ventional vidicons . In CCI (charge-coupled imagers), a 
scene incident on the MOS array surface is transformed into 
a distribution of minority charge packets under the gates. 

To read the image out, a sequence of voltage pulses are 
applied to the array which transfer the charge packets 
across the CCI and produce a video signal at the output. 

The device is smaller, much cheaper to build and simpler to 
operate . 
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Recently, Fairchild has reported a linear 5^0 element 
CCI array which has a 1,000 to 1 dynamic range and maintains 
good picture at 0.03 foot-candles.^ (Vidicons have about a 
100 to 1 dynamic range.) Such a performance is approaching 
the field requirement of low-light-level television (L^TV). 
Bell Telephone Laboratories, using 3 CCI chips of 106 x 128 
elements each, has demonstrated that CCI can be used as 
color TV cameras although much denser CCI are needed for 
practical applications (250 x 25 O elements for plcturephone , 
550 X 550 elements for commercial TV quality) [6]. There 
seems to be a good probability that CCI will replace most of 
the electron beam scanning devices in camera tube applica- 
tions. CCI can equally be used in other applications as 
solid state page readers, fascimlle sensors, information 
displays, computer terminal readouts, etc. 

(2) In anolog delay applications, CCD will also make 
an important impact replacing the present delay methods Ir 
signal processing using glass rods, quartz crystals or 
coaxial cables. In CCD, the signal is introduced and the.i 
transferred across the device by the clocking pulses. 

Delays from microseconds to milliseconds can be achieved by 
changing the clock frequency. For the first time, an 
anolog signal delay device can be fabricated by'IC process 
and Incorporated into other IC systems. The applications 



^Several other companies have achieved similar results. 
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are numerous in communications, radar, etc., such as the 
delay lines in TV sets, transversal delay line filters in 
telephone systems for equalizing phone lines, matching fil- 
ters for cooing a.nd decoding in radar systems, etc. 

(3) For digital memory, the Impact of CCD is not as 
immediate .mair.ly because of the highly developed magnetic 
core and semi conduct or memory technologies already in exist- 
ence and the other developing memory techniques such as mag- 
netic bubbles and optical holography. However, operational 
buffer mem.ory; consisting of two 480 bit CCD shift registers 
on a 110 rnil^ Si chip has been demonstrated [13]* For bulk 
memory, gross packing density of 3 million bits per square 
inch was predi'. ted . Cost of 0.01-0. let per bit is possible. 
It is concr i\alic that CCD may have stronger impact in 
memory than ;..n imaging and anolog delay. 

The interest of this thesis is in the application of 
CCD in Imaging, more specifically, in Infrared Imaging. The 
Interest in infrared imaging is quite broad indeed, such a~ 
the downward looking and forward looking infrared systems 
in military and security night vision missions, remote 
sensing in ecology and environment fields, nondestructive 
testing in industrial quality control and the very promising 
thermography in medical diagnostics. 

However, the feasibility of CCD concepts in IR had not 
been investigated at all at the start of this thesis (June 
1972). At it stands today, several approaches have been 
discussed. First, the monolithic MOS approach in narrow 
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energy gap semiconductors similar to the Si-CCD is the 
obvious extension to IR. Second, semiconductors of small 
energy gap are used as IR sensors, and the photo-excited 
carriers are transferred into the Si-CCD. Two subapproaches 
exist, one being the photoconductive sensor-Si~CCD package 
and the other being the IR seraiconductor-Si-CCD heterojunc- 
tion approach. 

This thesis considers the monolithic approach. It 
should be pointed out that another research program of this 
group, the development of polycrystalline narrow gap semi- 
conductor photoconductive IR detectors [l6][ll], can be 
applied to the photoconductor-Sl-CCD approach. In this 
monolithic approach, the first task is to examine if the 
surface layer of the new semiconductor in a MIS structure 
can be controlled by the gate voltage. Specifically, it is 
necessary to determine if: 

a. Accumulation, depletion and Inversion layers can 
be controlled by the application of proper gate voltages. 

b. Does the' surface layer stay long enough in its pc ’- 
tially depleted state before it becomes inversion layer? 

This is the so-called storage time which must be longer 
than the integration time in an imaging application. 

c. If the storage time is long enough, is the storage 
capacity adequate for handling the photon flux of a scene? 

Only after these questions receive affirmative answers, 
can we proceed to investigate the IR CCI applications. For 
this purpose, a group of I'V’-VI compound and alloy semicon- 
ductors are studied; PbTe, PbSe, PbSnTe and PbSnSe. Their 
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energy gaps are all less than 0.3 e^/ . wnich make them 
suitable for IR imaging beyond 4u, such as in the and 

8-12u atmospheric windows. MIS structures using AI 2 O 3 and 
SiOj as insulators and Hg probes as metal gates were fabri- 
cated and studied. 

The main purpose of this thesis is to investigate if 
depletion and in\ <;. : Ion layers exist at the interface of 
the insulator and semlccnductor in these MIS. If so, can 
they be controlled by the gate voltage? 

During the course of this research, MIS of two other 
semiconductors, Ge and Si, were also studied. The motiva- 
tion was to acquaint ourselves with the MIS research using 
MIS samples of two well known semiconductors. First, Ge- 
MIS using Al^G, as insulators were studied. The energy 
gap of Ge varies from O .67 ev at 300°K to 0.73 ev at 77°K 
which is suitable for I. 7 OU IR Imaging. It should be noted 
that Ge-MIS has not been studied before, making its study 
useful by its own right. 

Second, the effect of high-energy-electron bombardment 
on a n-channel Si depletion MOSFET was also studied. The 
electron energy was 65 .^ Mev. and the radiation was approxi- 
mately 5^10^^ electrons per square centimeter per second. 
Exposure to this radiation for five seconds is equivalent to 
exposure in the space environment around Jupiter for about 
ten hours. This study is also related to a research program 
at the Linear Electron Accelerator of the Naval Postgraduate 
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School studying the radiation effects on semiconductor 
devices during a space flight around Jupiter. 



19 



II. REVIEW OP CCD CHARACTERISTICS 



Since the charge-coupled concept is so new, a description 
of the device Is necessary. This section of the thesis will 
cover the concept, the structure, the method of operation 
and the "figures of merit" of the CCD. 

A. CONCEPT 

The basic concept of CCD's consists of storing minority 
carriers (or their absence) in a spatially defined depletion 
region (potential well) at the surface of a homogeneous 
semiconductor and to move this charge along the surface by 
moving the potential minimum. The potential minima are 
produced by applying a voltage to conducting electrodes 
(gates), formed on the surface of the Insulator covering 
the semiconductor, and driving the surface into depletion. 

A variety of functions can be performed by injecting charge' 
into the potential. well , transferring this charge along th 
surface of the semiconductor, and detecting the charge at 
some output location. 

B. STRUCTURE 

2 

The basic structure of a CCD is a linear array of MIS 



An MIS element is a metal-insulator-semiconductor 
structure whose main function is to produce a potential well 
in the semiconductor from a bias applied between the metal 
and semiconductor. MIS will be covered in detail in the 
next chapter of this thesis. 
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elements as shown in Figure 2-1 [1]. Figure 2-1 also de- 
picts the storage and transfer schemes in a CCD. In Fig. 
2-la, every third electrode is biased more strongly than the 
other two and charge can be stored in the potential well 
under gate 1. In part b., the transfer sequence, the bias 
on gate 2 is made greater and as a result, if the gate 
spacing is close enough to allow adequate coupling, the 
charge transfers into the region under gate 2. Part d. in 
Fig. 2-1 indicates the clock pulse scheme used to transfer 
the charge. In this case the device is a 3 phase ( 34 >) CCD 
since three different voltage levels are required to move 
the charge. Other type schemes are described below. 

C . MODUS OPERANDI 

The operation of the CCD is based on the minority 
carrier charge formation, storage, and transfer in an inver- 
sion layer at the semiconductor surface. The composite 
element of a CCD in the case of the 34 > device are 3-MIS 
basic elements. There are, however, other clocking scheme ; 
that are used. Note also that the charge transfer can be 
in any direction dependent on the nature of the voltage 
pulses. Inherent disadvanteges in the 34) device are the 
criticalness of the interelectrode spacing, cross unders 
required in device fabrication, and the gap areas between 
electrodes which can assume potentials that adversely 
affect the transfer process. 
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Figure 2-1. Charge-Coupled Device Basic Structure and 
Transfer Process: (a) Cross section of the 
charge-coupled structure, (b) surface potential 
profile at the beginning of the transfer cycle, 
(c) surface potential profiles during the trans- 
fer cycle, and (d) three-phase clock waveforms. 
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Figure 2-2 shows the construction and operation of the 
2(f)-CCD. This type CCD not only Is simpler to construct 
than the 3< but also has the advantage of minimizing the 
separation between electrodes. The direction of charge 
transfer is determined by the asymmetry of the potential 
wells under tv;o adjacent electrodes; the smaller potential 
well beneath the 'thicker oxide causes the packet to flow to 
the right only. The asymmetric potential wells necessary 
to produce directionality can also be obtained v;ith an ion 
implanted barrier as shown by the shaded region in 
Figure 2-2. 

Other CCD structures include -CCD devices and a 
Uniphase CCD just recently proposed [9] in which a silicon 
nitride layer with a built-in charge allows one-directional 
charge-packet flow. 

The input and output schemes of the CCD are depicted in 
Fig. 2-3. Optical generation, the scheme most important i". 
this thesis project, is listed in part c. of Fig. 2-3. 

Note that in the figure the radiation is introduced in tht 
back of the device, through the semiconductor. Other de- 
vices pass the radiation through windows in the front. For 
example, Melen [17 J uses polysilicon gates between aluminum 
gates through which the light image passes. Other devices 
coming into being are CCD with depleted-buried channels in 
which the charge packet flows. These devices will gain 
rapidly in popularity because of the elimination of surface 
states which plague surface CCD's. 
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Figure 2-2. Construv, -fi .. id Operation of an Ion Implantation 2<j)-CCD 
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Figure 2 - 3 ‘ Input and Output Schemes of a CCD 



D. FIGURES OF MERIT 



Several figures of merit have been defined: quantity of 
stored charge, efficiency of the transfer, equilibrium time 
after deep depletion (storage time) and charge conservation. 

1 , Charge Storage Capability 

This is the ability of the CCD to store charge. It 
depends on the following factors: energy gap, doping and 
dielectric constant of the semiconductor, the dielectric 
constant and the thickness of Insulator, and the interface 
states of the MIS structures. A measure of charge storage 
capacity is the capacitance of the device which can be 
obtained from the C-V relationship. The stored charge is 
divided between the insulator and the depletion capacitance 
in the semiconductor. It is, of course, beneficial to have 
most of the charge in the semiconductor, however, there is 
a maximum charge at which the surface potential of a cell 
has the same value as its adjacent cells and the device can 
no longer function. A useful plot in estimating maximum 
charge ability is a plot of surface potential vs applied 
gate voltage. 

2 . Transfer Efficiency 

One of the most important figures of merit is the 
amount of charge transferred or the transfer efficiency, n. 
Often the efficiency, or lack of it, is described in terms 
of the percentage of charge remaining after a transfer. The 
efficiency can be discussed in terms of element efficiency 
per transfer or device efficiency and values up to 99-995^ 
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are common [26] . The importance of •:his figure of merit is 
amply demonstrated in the various reports on the worth of 
CCD's I'/here the transfer efficiency is the prime interest 
[3, 7, 17 , 27 ] . A measure of the difference between the 
amount of charge that exits a CCD to the charge input by any 
means is given by n. The cause of this loss is of major 
concern . 

One of the major sources of charge loss are the 
Interface states. These states are energy levels within 
the energy gap of the material that fill up rapidly but 
empty more slowly ana as a result some travel with the 
charge packet from v/hence they came and others join up with 
a following packet. The study of Interface state density in 
Pbi_xSn^Te (Se ) alloy semiconductor MIS samples is of major 
interest in this thesis project. Two ways to eliminate 
this detrimental effect is to fill the Interface states 
with background charge (the so-called fat zero approach) o ■ 
move the interface deeper into the semiconductor (the so- 
called buried channel approach) . 

3 . Transfer Time, t ^ 

The transfer time, defined by the time necessary for 
charge transfer to an adjacent cell, is due primarily to 
diffusion but is aided significantly by parallel’ interelec- 
trode electric fields existing in the device. The transfer 
time should be shorter than the thermal diffusion time for 
minority carriers. 
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• storage Time, t„ 

The storage time is defined as the time for a pulsed 
CCD element to return to the steady state conditions. It 
is a function of device capacitance, applied voltage, and 
the dark current of the device. 

Storage time is an Important parameter since it 

determines the lower limit of the operating clock frequency 

of the device. Storage time must necessarily be longer than 

the period of the clock frequency or the device will not 

function properly. This thesis project will also consider 

storage time which can be measured in principle. Typical 

values are t >0.5 sec for Si-MOS-CCD and for Imaging applica- 
s 

tlons t must be much greater than the optical integration 
s 

time. Storage times are affected by surface states whose 
energies lie within the energy gap of the semiconductor. 
Surface states consist of layer states characteristic of 
the oxide and interface states. Of the various types of 
layer states, interface states (of lower density) are more 
significant since they have a higher minority carrier cap- 
ture cross section. 

5 . Charge Conservation 

Charge can be lost or gained in the Interelectrode 
gap region. If the region between electrodes is' allowed 
to enter accumulation, and there are minority carrier traps 
present, recombination and trapping will occur and minority 
carriers are pumped into the substate and lost. To prevent 
this loss, one must ensure that the interelectrode gap 
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regions are depleted. The reverse may occur if minority 
carriers from the substrate move into the charge packet at 
the edge regions of the gate electrodes. This effect is 
small if channel-stop regions arc- usee. 

E. NON-INSULATOR CCD 

The surface states are a real concern in CCD’s and a 
possible way to eliminate this adverse effect is to provide 
a buried channel in which the minority carrier charge can 
form and move. The buried channel essentially eliminates 
the surface state trapping effect. A uniquely different 
buried-channel CCD which has no oxide interfaces has recent- 
ly been proposed [22] . This section is presented separately 
for reasons vjhich will be evident later. These reasons are 
the result of evaporating a form of alumina which had 
characteristics of a conductor and not an Insulator. 

The two proposed CCD’s are the double-j unction CCD 
(DJCCD) and the Schottky-barrler CCD (SBCCD); both CCD’s are 
void of an insulation layer, characteristic of CCD’s. Thf 
DJCCD is a p"*"-n-p structure for n-channel CCD and the 
potential well is formed between two reverse biased pn Junc- 
tions. Figure 2-4 gives the conduction band profile of the 
device. Curve A in the figure applies for OV on the gate, 

B for +2V and C for -2V. The potential well is clearly 
formed in the n-type region and, as can be seen from the 
figure, more that 90 % of the potential applied to the device 
goes into providing the depth of the potential well. This 
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Figure 2-4. Theoretical Conduction-Band Profile for a 
Double-Junction CCD 
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effect allows for excellent efficiencies as fai’ as charge 
capacity is concerned. 

The SBCCD is basically similar to the DJCCD. However, 
the p"*" region is replaced by a metal layer. Its character- 
tics are similar to the DJCCD except that leakage current 
is larger. 

These devices are very promising because of various 
reasons, the most important reason is the elimination of 
the interface states with the removal of the insulator 
layer. This reason is a distinct advantage to the labora- 
tory that can produce good semiconductors but bad insulators. 
A second real advantage for the DJCCD is the fact that the 
clock driving voltage can be positive, negative or sinusoidal 
since the depletion region is produced by varying the dc 
bias on the gate and not on the clock pulse. A third ad- 
vantage is that these devices are not so vulnerable to radia- 
tion as standard MIS CCD's are since the radiation affects 
the insulator greatly and the insulator has been eliminatec . 

The feasibility of producing such devices in the Naval 
Postgraduate Solid State Laboratory is somewhat promising 
since good semiconductor deposition, on Pb-Sn-Te(Se) alloys 
has progressed to a fine art. On completion of the multi- 
source deposition scheme, the laboratory will have the 
capability to produce such devices on a daily basis. Hetero- 
junctions have already been produced for application to 
heterojunction lasers [10] . 
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III. REVIEW OF MIS CHARACTERISTICS 



The basic characteristics of the CCD are deterir.ined 

3 

by studying the MIS. Three major parameters determined 
are the surface potential, 4>g, -a function of applied 
voltage, the storage time, and the surface state density. 

The surface state density affect: the transfer efficiency, 

the storage time puts a lov;er limit on the clock frequency 
of the CCD and dg is an indication of the depletion width 
which is related to the amount of charge stored. This chap- 
ter first reviews the basic properties of the MIS and the 
methods of their investigations sue' as the capacitance- 
voltage (C-V) curves, the conducts: ce-voltage (G-V) curves, 
the capacitance-time (C-T) curves and the derivatives of the 
C-V and G-V curves. This discussion pertains both to the 
ideal MIS which contains no interface states (Ng^) , no metal- 
to-semiconductor work function difference ('t’j/jg)j negligible 
bulk semiconductor and insulator resistance, and the non- 
ideal case where these parameters are considered. 



A. DESCRIPTION OF AN IDEAL MIS 

In this section the basic MIS structures are discussed. 
It is important to note that no nev; theory was developed in 



3 

The MIS structure is the basic cell of a CCD whereas 
the basic element of a CCD is a set of 2, 3 or ^ cells (MIS) 
depending on the clock pulse scheme used. 
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this thesis but, rather, MIS structures of lead-tin alloys 
were analyzed following the existing theory developed for 
silicon MIS devices. 

Figure 3-1 provides a diagram of the energy bands and 
charge distribution in an ideal MIS for the accumulation, 
depletion and inversion conditions in a p-type semiconductor. 
The development here will deal mostly with p-type devices 
since p-type lead-tin alloy semiconductors were in abundant 
supply to be tested whereas the supply of n-type devices 
was limited.^ 

5 

In a p-type MIS, when a negative bias (positive bias 
for a n-type MIS) is applied to the gate, the majority 
carriers (holes) are attracted to the interface between the 
semiconductor and Insulator by the attractive forces of the 
net negative charge on the gate. Accumulation is depicted 
in part a. of Fig. 3-1 with the bands bending upward. With 
no bias applied, a flat-band condition exists in which the 
energy bands are flat. When a small positive bias is 
applied to the gate, majority carriers are pushed away from 
the surface and a depletion region is formed. This is 
shown in part b. of Fig. 3-1 and the bands are bent 



4 

The prepared n-type sample concentrations were in the 
mid 10^® range and showed no C-V variation. 

5 

The bias on a MIS will always be with reference to the 
substrate . 
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Figure 3~1* Energy Bands and Charge Distribution in Ideal MIS for 
Various Bias Conditions 



downward. If the positive bias is increased further, 
minority carriers (electrons) will be attracted to the 
interface and form an inversion layer. As the positive 
bias becomes even larger, more and more minority carriers 
will exist in the inversion layer. In the latter two cases 
the energy bands are bent down further as shown in part c. 
of Figure 3-1. 

It is important to note that the above description, in 
the depletion and inversion region, applies only if the bias 
change is slow enough for minority carrier formation at the 
interface. This condition is the dc or low-frequency case 
in which it is assumed that the minority carriers are 
generated and moved to the surface and that the majority 
carriers (holes) formed in the rest of the depletion re- 
gion are removed. The minority carrier charge formation 
may not always occur. Vfnether it does or not is of major 
Importance to CCD applications since the rate at which the 
Inversion layer appears determines the storage time. In 
the development v;hich follows it is assumed that the dc 

bias is varied slowly enough for the minority carriers to 
7 

follow . 



The Intrinsic condition occurs when an equal number of 
minority and majority carriers exist at the surface. The 
mid-gap energy level at the surface, in this case, is equal 
to the Fermi band, and the corresponding bias voltage is 
referred to as the threshold voltage. 

7 

This concern does not apply to the majority carriers 
whose relaxation time constant is very short, t = e/a. 
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An expanded view of the energy band diagram of a p-type 
MIS under inversion condition is shown in Pig. 3-2. The 
parameters used in Figux’e 3-2 are: 



Ey = valence band level 
= conductor band level 
E^ = intrinsic band level 
Ep = Fermi level in semiconductor 
Ep^= Fermi level in metal 
4>sc“ T^ork function in semiconductor 
~ function in metal 

4> = potential in semiconductor 

= potential in semiconductor bulk 

4> = potential at semiconductor surface 

In the ideal case, the space charge within the semiconductor 

is given by 



= / p(X)dX (1) 

^ X=0 

where p(X) = q (p-n + - N^) , X = 0 is at the interface, 

p and n are the carrier concentrations and is the 

Impurity concentration. For nondegenerate semiconductors 
in thermal equilibrium, Eq. 1 can be solved by applying 
Poisson's equation and Gauss’ law in the space charge region 
[4] . The solution is 

u f 

^s ^ jiT^ qnj_Lj3 1^2 [cosh (u^-u^) -cosh u^^ 

+ Ug slnh (2) 
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Figure 3-2. Expanded View of MIS Energy Band Diagram Under 
Inversion Conditions 
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where is the Debye length given by 



■'D 



f k r V- 

q 2qn^ 



( 3 ) 



and in kT/q units, is the potential at the surface 

(in the bulk) of the semiconductor. Q^, which is the inver- 
sion layer contribution to Q , is given by 

s 



-u 

s T 

v-\ 1 






f exp (u - u^) du 


FIT " “ 




1 


1 2[cosh(u-u^ )-cosh u^+u sinh u^]'^ 



(4) 



These quantities are important since from Qg the space 
charge capacitance can be determined and from the charge 
in the inversion layer can be determined. The latter quan- 
tity is instrumental in describing the charge storage capa- 
city of the device. 



B. CAPACITANCE-VOLTAGE CHARACTERISTICS OF AN IDEAL MIS 
One of the best known and widely used tools to study 
the nature of the MIS is the capacitance-voltage (C-V) 
curves. Shown in Fig. 3-3 is an ideal C-V curve for a 
p-type MIS. Notice from the figure that the vertical scale 
contains capacitance normalized to the insulator capacitance 

(C^). 

1 . Description of Typical C-V Curves 

The information in a C-V plot is contained in the 
depletion and inversion regions of the capacitance plots. 

In the accumulation region the maximum capacitance the MIS 
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Figure 3-3. Ideal C-V Curves for p-Type MIS 



device can have is the insulator capacitance. The insulator 
capacitance shields any detectable change in surface 
capacitance . 

I 

The capacitance in a MIS structure can be divided 
into 5 types. These are the accumulation, depletion, low- 
frequency inversion, high-frequency inversion, and high-fre- 
quency deep depletion capacitances (depicted in Pig. 3-3)* 

The first three types of capacitance were described 
in section A above. In the high-frequency inversion case the 
ac signal frequency is so high that the change of inversion 
charge (increase and decrease of minority carriers due to 
ac signal) can not be supplied by the generation and recom- 
bination rates in the semiconductor and the capacitance lev- 
els off. This condition is really a partial inversion case 
where the minority carriers act like fixed charges and the 
semiconductor capacitance is effectively given by = 

£ /X. and exists for surface potentials greater than 
about two times the potential in the bulk. 

In the high-frequency deep depletion case the dc 
bias sweep and ac signal are both rapid enough not to allow 
the inversion charge to form. As a result, the capacitance 
value decreases even further from the Inversion case and 
goes to zero, in theory. 

2 . Theoretical C-V Equations 

In the low-frequency inversion case the total capa- 
citance of the device is given by the series combination of 
the insulator capacitance and the space charge capacitance. 





1. J 
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which when normalized becomes: 
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1 + C./C 3 
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Since C_ = - -r-r^ 
S o<j) 



( 5 ) 



, differentiation of Eq. 2 



for 



Q above and substitution into Eq. 5 will yield the expres- 
s 

Sion for the ideal low-frequency C-V value. 

For the high-frequency inversion case, an effective 
depletion depth must be considered since the minority 

carriers do not follow the ac signal but are affected by it. 



The space charge capacitance per unit area is defined as 
K. 



-s ^0 

C„ = Y where is X^-_ 

S 6X1 0il 



U 



u 



1 - e 
F(u,u. ) 



du 



( 6 ) 



and P(u,Uj^) is the term enclosed in brackets in Eq. 2. The 
high frequency capacitance then becomes Eq. 5 with Eq. 6 
substituted for Cg. 

In high-frequency deep depletion case, the C-V rela- 
tionship is like that of a reverse biased p-n junction. 

This is really a transition case in which the MIS is not in 
equilibrium and can be caused by a leaky oxide or if the dc 
bias is rapidly swept from accumulation into inversion. 

The capacitance is again given by Eq. 5 with a different 

value for X since in this case there is no inversion 
eff 

layer and Q = q(N_ - N )X . The capacitance can be found 
S l) a 0 X X 

to be [4] 
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The deep depletion capacitance value Is important in des- 
cribing the storage time of the MIS. 

C. CAPACITANCE-TIME CHARACTERISTICS OF MIS 

When a MIS structure is pulsed from accumulation into 
deep depletion, the semiconductor surface potential is large 
at the onset of depletion but relaxes to the quasi-equilib- 
rium Inversion value in a certain characteristic time, T. 
This time is dependent on the rate of minority carrier 
generation at the interface and the rate at which the major- 
ity carriers drift to the bulk region to neutralize the 
ionized acceptor atoms. This C-t response characterizes the 
diffusion length, the minority carrier lifetime, and the 
surface generation velocity of the device. Other Important 
parameters derived from the C-t response are the storage 
time, T, and a measure of the maximum amount of charge that 
can be stored in a MIS under pulsed conditions. 

The C-t response of a MIS exists because the minority 
carriers can not follow the ac signal (or the pulsed dc 
bias). Figure shows the relationship between the C-t 

and C-V curves. After the MIS capacitor is pulsed into 
deep depletion, it relaxes to the quasi-equilibrium high- 
frequency inversion in time T. This relaxation is due to 
the appearance of minority carriers at the Interface by five 
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Figure 3“^ • Accuirmlation to Inversion Pulsed Response of an MIS 



possible generation mechanisms [21]; thermal generation in 
the depletion region, thermal generation in the semiconduc- 
tor bulk region, thermal generation at the surface of the 

bulk, and external excitation such as photons and energetic 

\ 

electrons. The consequence of these mechanisms is discussed 
below. 

1 . Storage Time 

The time, T, is a measure of MIS storage time which 
determines the maximum amount of time that a MIS, once 
pulsed, remains in deep depletion. In MIS-CCD application^ 
this time determines or sets the lower limit on the clock 
frequency . 

In the absence of an external source of minority 
carriers, storage time is controlled by the rate of thermal 
generation of minority carriers. In semiconductors, like 
the indirect gap silicon, this is a slow process and may 
take up to 200 seconds. Long storage times are favorable 
in CCD applications since it allows fewer unwanted non-sig- 
nal carriers to enter into inversion. If the charge is 
dominated by the thermal generation in the depletion layer, 
this storage time is approximated by T = 2tgN/n^ where 
is the minority carrier lifetime, N is the impurity concen- 
tration, and n. is the intrinsic carrier concentration. For 
practical CCD operation, the clock period of transfer should 
be at least one hundred times shorter than the storage time. 



2 . Storage Capacity 



An Important consideration for the CCD operation is 
the maximum amount of minority charge that can be stored in 
each MIS element. This maximum, charge can be thought of as 
the amount of charge v/hich when placed beneath the gate in 
the depletion region area of a MIS element in a CCD will 
raise the potential to that of its adjacent neighbor (to 
that of the surrounding area) i.e., the surface potential of 
the transfer element becomes that of the storage element. 

In this thesis, the storage capacity is estimated 
by the following procedure. First, the theoretical C-V 
high-frequency inversion and high-frequency deep depletion 
capacitance-voltage curves are computed. Q storage is esti- 
mated as Q , = (C. - C ) V. where V. is the gate 

store inv DD inv inv ^ 

g 

voltage in the inversion region. 



D. OTHER RELATED MIS CHARACTERISTICS 

Other methods used to investigate a MIS structure are 
the conductance-voltage (G-V) curves and the derivatives r r 
both the C-V and G-V curves. The G-V curves are not as 
popular as C-V curves in describing a MIS because not as 
much information can be obtained from these curves as can 
be obtained from the C-V curves. For instance, the dC/dV 
curve’s peak is a good indication of the flat -band voltage 
and G-V curves allow the determination of interface state 



g 

°An exam.ple of the output obtained from the computer 
program output appears in Chapter VIII. 
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density. The Information obtainable from the C-V curve 
are carrier concentration, surface potential versus applied 
voltage, and interface states among others. The dC/dV curve 
Is similar in appearance to the G-V curve in that both have 
peaks in the depletion region and are flat in accumulation 
and inversion. The experim.ental determination of these 
curves was not attempted in this thesis. 

E. C-V CHARACTERISTICS IN REAL MIS 

The shape of *ohe ideal C-V curve, shown in section B, 
above, pertains to a structure with no insulator charge, no 
work-function difference between metal and semiconductor, 
no interface states, and no surface charge. Flat-band condi- 
tion (l.e., where the surface potential goes to zero) is 
also zero in this case. In this section the characteristics 
of a real MIS are discussed. 

1 . Physical Changes 

In Fig. 3-5 the changes in the real insulator which 
made it different from the ideal case are shown. A first 
change is the insulator charges which could include the 
fixed Insulator charges, ionized traps, mobile ions and in- 
terface states. Second is the interface states in the 
'semiconductor. Third is the work function difference. 

The insulator charge has the effect of producing a 
change in the potential of the semiconductor. This effect 
was deliberately introduced in the oxide Insulator of sili- 
con nitride in the unlphase CCD described in Chapter II, to 
produce potential wells in the semiconductor for directional 
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Classification of Charges and States in a Non 
Ideal MIS 



charge packet flow. The Interface states within the for- 
bidden gap of the Insulator and semiconductor ("traps") 
are occupied and emptied at different rates. Fast Interface 
states react quickly and are usually found at the semicon- 
ductor-insulator interface whereas the slow states react 
slowly and are found at the metal-insulator Interface. 

2 . Device Characteristic Changes 

These changes refer to changes in the device charac- 
teristics such as the C-V, G-V, C-t and their derivative 
curves. In addition to the above mentioned causes of non- 
ideal behavior, the following factors also effect these 
curves. These are high contact resistance and high bulk 
generation lifetime. 

When the only departure from ideal is the work 

function difference , the C-V curve is simoly shifted 

ms * 

parallel to the left or right and the amount of shift is 

9 

equal to this difference. liiHien the departure includes 
interface states that are Independent of the swept dc bla , 
then the curve is also a simple parallel shift on the volcage 
axis. The amount of interface states is a measure of the 
charge occupancy of the states. The voltage shift including 
these two effects is AV^ = - Q-./C.. 

(jr inb bo 3- 

When the interface states become dependent on the 
applied dc bias, whose density varies 'within the energy gap. 



^The shift is to the left for negative (f)„_ 
right for positive <i>jvrc • ^ 



and to the 
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the curve Is shifted in a non-parallel fashion along the 
voltage axis but the ratio of maximum to minimum capacitance 
is unchanged from the ideal case. In the most extreme case 
where the surface states do not follow both the ac signal 
and dc bias, the maximum amount of deviation from the ideal 
case exists. Figure 3-6 provides a family of curves showing 
the effect of the non-ideal parameters described above. 

These curves can be considered a ’’finger print" of non-ideal 
curves and as such can be used to rapidly evaluate the 
quality of a MIS structure after an experimental curve is 
obtained. Note from the figure that part a. is simply the 
deep depletion case described in section B. 

Also displayed in Fig. 3-6 is the hysteresls^*^ 
effect present in some MIS samples. The amount of hysteresis 
in C-V curves is determined by the amount of current passing 
into and through the insulating film and by the amount of 
charge that is trapped in the insulator. Device hysteres: 3 
does not necessarily make it useless. It is feasible to 
obtain a greater quantity of charge storage in devices di .- 
playing hysteresis. 

Another non-ideal characteristic is drift or the 
Inability to trace repetitively similar characteristic 
curves. 



Hysteresis is the inability of a particular MTS sample 
to characterize similar C-V curves when the device is biased 
from acc'amulatlon to Inversion or from inversion to accumu- 
lation. Hysteresis also applies to G-V, I-V and the deriva- 
tives of C-V and G-V curves as well. 
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Figure 3-6. "Finger Prints" of Non-Ideal MIS 
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3 . Theoretical Changes 

Two of the major non-ideal effects 
in the potential distribution in the MIS. 
fied in the Eqs. 8 and 9 • 



= h " '"s " *MS 



are the changes 
They are signi- 

( 8 ) 



-Q 



M 



= ^1 






Qc 



(9 ) 



with the inclusion of the work function difference the 

insulator charge, and interfacial charge. Plat -band voltage 
is also changed and does not occur at V„ = 0 but is defined 
“ *^ss' ’ where is the overall effect by 

the added charge. These changes can also be represented in 
the admittance of the MIS as shown in the equivalent circuit. 
Fig. 3-7. It can be obtained for experimental C and G 
measurements [28]. Their expressions are long and are not 
Included in this thesis. The combined device capacitance 
in low-frequency and high-frequency cases become 



C = 



"i 



(C< 



Coc.) 

oo 



^s ^ss 



(low frequency) 



(lO ) 



and 



C = 



C. 

1 



C. 

1 



+ c. 



(high frequency) 



( 11 ) 



These two equations are used in the computer program to 
analyze the experimental data. Equation 11 was used above 
to calculate both the theoretical high and low-frequency 
capacitances with different values of for each respectively. 
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Figure 3-7* MIS Equivalent Circuit 
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IV. MEASUREMENT TECHNIQUES ON MIS 



The basic schematics of the equipment used and the 
problems associated v;ith the experimental procedure are 
described in this chapter. Three measurements were used: 
high-frequency C-V, low-frequency C-V and high.-frequency 

C-t . 

A. HIGH-FREQUENCY C-V 
1 . Basic Schematic 

The basic schematic for obtaining experimental high- 
frequency plots is shown in Fig. 4-1. The Boonton Model 72A 
capacitance meter provides a direct scale readout on X-Y 
recorder output for full-scale capacitance ranges of 1, 3, 
10, 30, 100, 300, 1000, and 3000 picofarads. It requires 
a minimum of setup time. The bias voltage which ranges from 
i200, i400, or ±600Vdc is applied to the back terminals 
and automatically appears on the front test coaxial input 
connectors along with a 15 MV-lMHz ac signal. The meter's 
phase sensitive detection system allows accurate capacitance 
measurement of devices with a Q as low as one. A dc output 
voltage which is proportional to the meter capacitance read- 
ing is available at output Jacks for the input to an X-Y 



The bias signal appears on the low-coaxial input con- 
nector if it is applied to the lo back Jack input and vice 
versa. 
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Block Diagram for High-Frequency C-V , and 
Measurements 



Figure 4-1. 



recorder in terms of 1 or 3 maximum output voltage for the 1 
and 3 series scale settings respectively. 

2 . Errors Due to Leaky Insulators 

During the first phase of experimentation, difficul- 
ty was encountered with the measured capacitance of the 
lead-tin samples. The capacitance meter gave some readings 
on one scale but not at all on other scales. For example, 
the capacitance was 28 picofarads on the 30-p scale but over- 
loaded the scale on the 100 p-scale. The cause of this 
difficulty was not due to the Boonton capacitance meter 
since it was pre- and post-tested for accuracy. 

The first cause investigated was the possibility 
thdt the mechanical probe punched through the metal gate 
and the insulator. A second cause investigated was the pos- 
sibility that the high electric field existing at the tip 
of the probe produced dielectric breakdown. To avoid these 
problems, soft indium solder on the probe tip was used. 

This prevented the high E field at the probe tip by spread- 
ing out the probe tip over a larger area. It did not elimi- 
nate the problem. When the mechanical probe tip was 
replaced by a small dot of liquid Hg, the problem among 
different scales was eliminated in som.e cases. 

Further investigation showed that the dc leakage 
resistance of the Insulator along the Insulator surface 
and across it was approximately 1 Kfi . This check provided 
a rapid evaluation of the sample as to whether it would 
produce useful C-V curves. Only those samples whose 



55 



leakage resistance was greatei"* than about 2YSI gave useful 
curves . 

In order to verify the effect produced by a leaky 
insulator, a shunt resistance box v/as placed across the MIS 
samples which were known to have a high value of dc resist- 
ance. Figure 4-2 and 4-3 show the affect produced for 
various values of shunt resistance. It is clear from the 
figures that shunt resistance (leaky Insulators) produces 
a drastic change in the C-V curves of a MIS sample. In 
Fig. 4-2 the sample was a silicon MOS ’with insulator resist- 
ance of greater than 1 MQ. In Fig. 4-3 the same effect 
as Fig. 4-2 was produced on a Ge~MIS which had an insulator 
resistance of 10 It suggested tha.t the C-V curves for 

a MIS with leaky Insulators less than 1 resistance do 
not produce meaningful results. For leakage resistance 
between IKJ^ and lOKJ^, C-V curves are probably distorted. 
However, the changes in C-V depletion and inversion is 
still clearly evident. 

B. C-t CURVES 

1 . Basic Schematic 

C-t measurements were carried out using the same 
setup as for the C-V measurements. A square wave from the 
Wavetek signal generator was used. The voltage variation 
necessary to switch from accumulation to Inversion was ob- 
tained from a check of the C-V curve which was accomplished 
prior to the C-t measuremient . For those MIS whose storage 
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Figure 4-2. Leaky Insulator Effect on n-Type Si-MOS with Insulator 
Resistance Greater than iMfi 
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Figure ^-3- Leaky Insulator Effect on p-Type Ge-MIS with Insulator 
Resistance of liKQ 



times are too short to be recorded on the X-Y recorder, a 
Tektronix Model 422 oscilloscope was used to observe the 
C-t response. 

2. Errors Associated with C-t Measurements 

C-t curves are unmeasurable when the insulator is 

leaky because of the available supply of minority carriers. 

For this reason, m.ost lead-tin MIS samples displayed no C-t 

response. Moreover, C-t curves are observable only if a MIS 

12 

displays a high-frequency deep depletion capacitance. If 
no characteristic deep depletion capacitance is observed, 
the semiconductor is already in the quasi-equalibrium inver- 
sion and can not relax any further. 

C. LOW-FREQUENCY C-V CURVES 

Because the Boonton test frequency is set at 1 MHz, a 
different experimental setup is required for low-frequency 
C-V curves. The C-V experimental test setup is depicted in 
Fig. 4-4. 

The Princeton Applied Research (PAR) Model 124 Lock-In- 
Amplifier (LIA) is the critical part in the experimental 
setup. A thorough understanding of how it functions is 
necessary before it can be used efficiently in obtaining 
meaningful experimental C-V curves. This section will 
consist only of a brief explanation of the LIA. 



12 

High-frequency deep depletion capacitance depends on 
the sample type and also on the conditions under which the 
sample is measured, l.e., temperature and illumination 
conditions . 
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Figure 4-it. Block Diagram for Lov;-Frequency and G-V 
Measurements 



60 



I 



I 



Considerable Information on its operation can be obtained 
from the various PAR application notes. [5, 25]. 

The LIA's forte is its ability to "lock-in" the signal 
frequency to the center frequency of the pass bands of the 
measuring circuits and tuned amplifiers in the device. This 
"lock-in" feature can be obtained for an external oscillator 
used as a reference or it m.ay lock-in internally, depicted 
in Pig. 4-4, to an internally generated sine wave signal. 

The signal detection ability extends down to the nano- 
volt range but at the high-sensitivity scale required for 
these nanovolt measurements, the noise problem is consider- 
able. Fortunately the required sensitivity scale for C-V 
curves was high enough to eliminate the high sensitivity 
noise problems. 

Signal detection in the LIA is accomplished through the 
use of a phase sensitive detector. The phase detector pro- 
duces an output that is proportioned to the signal input 
amplitude, /V /, and the phase difference between the refer- 
ence and signal input, <})o-4) . The magnitude of the refer- 

ence phase is not considered in the measurement. 

Figure 4-5 shows the equivalent circuit which describes 
the low-frequency phase sensitivity measurements, V, V^, 
and Vg are the phasor quantities of the voltages at the 
various positions shown. When these four voltages are com- 
bined with the Thevenin equivalent voltage of the oscillator 
output as seen from the gate of the MIS, the follovring 
equation results 
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Figure M-5- Equivalent Circuit for Low-Frequency Phase 
Sensitivity Measurements 
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G^((G-!-G^)cos (?-coC sin b) 



AV 



Vo = 



th 



(G+G^)2+to2(C+C^)2 



+ (o)(C+C^)cos B + G sin B)| 



( 12 ) 



where Y = G + juC is the admittance oi the substrate of the 
MIS to ground including the amplifier. 

the admittance of the MIS and B ft' tf-'^ ootal phase angle 
given by, B = (})^ - 4)^ - 4)^^ * 

Equation 12 simplifies to the follov;ing equation when 
G>>G^ and C>>C^^^ 



Vo = 



AV 



th 



(13) 



Thus, the output is a measure of capacitance if the two, 
restrictions on and G^ are m.et. A similar result obtains 
for conductance. 

The accuracy of Eq. 13 depends on the proper adjustment 
l4 

of 3 . If 3 deviates from the correct value by an amount 
A3 then the capacitance deviates from its true capacitance 



A^ ^ 

Cx “ toC^ 



A3 



(lit) 



where the quantity “ *^X quantity of the MIS. 

For a capacitance of 1000 pf, a resistance of 10 and a 



^^Thls relationship exists if G is on the order of 10 
and C is about lOOOp . [28]. 

^^3 is controlled by the LIA phase knob through and 
by the amplifier sensitivity setting through (p ^ . 
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frequency of 1 KHz, Eq. 1^1 becomes 

AC . 

= 100 A3 . (15) 

Thus the setting of 3 is critical in obtaining error-free 
capacitance curves. On the other hand, it can be shown 
that ^ -010A3 and apparently the conductance can be 

measured with little error for relatively large phase errors. 

Two types of preamplifiers are useable with the LIA for 
C-V measurements. The Model ll6 differential preamplifier 
requires a measuring resistor across its input to ground. 

The value of this measuring resistor depends on the conduct- 
ance of the MIS under test. For MIS conductance of 10 
the measuring resistor should have a value of approximately 
lOOfi. For larger conductances, the measuring resistance 
also increases. A measuring resistance value of was 

required on the Ge-MIS. The Model 18^ photometric preampli- 
fier is the better of the two amplifiers for three reasons. 
These reasons are: 1) an input signal monitor is provided, 

2) it is a more sensitive current sensing amplifier which 
requires no measuring resistor, and 3) absolute capacitance 
values can be obtained by calibrating the LIA with standard 
capacitors. Since the l84 has a virtual ground input im- 
pedance, it can not be used with devices that show high 
leakage . 

The LIA can also be used to give high-frequency C-V 
since it is sensitive to ac signal frequencies up to 210 KHz 
and this frequency is the high-frequency capacitance case 
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on many MIS (Si-MOS, for example). It can also give con- 
ductajice measurements by simply changing the reference 
phase angle by 90°. The derivatives of the C-V and G-V 
curves are also obtainable with slight; modification of the 
test circuit but were not attempted in dills thesis. 
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V. ANALYSIS OF C-V DATA 



Two computer programs were used to analyze the experi- 
mental C-V data in this research. They were originally 
developed for silicon MOS. It is assumed that they can be 
applied to germanium and the IV-VI alloy semiconductor MIS. 

A third program was developed to calculate the high-frequen- 
cy inversion capacitance minimum for different values of 

semiconductor dielectric constant and different type 
15 

Insulators. 

A. COMPUTER PROGRAM TO ANALYZE THE HIGH- AND LOW-FREQUENCY 
C-V DATA 

This program calculates and plots the following types 
of information from experimental high- and low-frequency 
C-V data: the impurity concentration at the semiconductor 
surface and within the bulk, plots of the variation of sur- 
face potential with respect to the gate voltage (V„-V ), the 
inputed C-V data (C-V), the inverse normalized space charge 
capacitance squared versus surface potential ((C . /C^ ) ^ -U_) 

the interface state density versus the semiconductor poten- 
tial (N 33 - 4 ')j and the effective interface charge versus the 
potential drop across the semiconductor surface (Qg^f^-Vg). 

15 

"^The dielectric constant of the IV-VI alloy semiconduct 
or samples is not known exactly. Doshler, in his theoreti- 
cal calculation of C-V, used 400 for this constant [8]. 
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This program was adapted from a program developed by the 
University of New Mexico for analysis of Si-MOS devices for 
radiation vulnerability studies [28], 

1. Semiconductor Surface and Applied Potential 

This section and the following tvjo sections contain 
only those equations necessary to allov; an understanding of 
the computer program and, as such, the equations will not 
be developed in detail. Consider the differential capaci- 
tance of the MIS 



C = 




( 16 ) 



where Qm is the charge on the gate and V is the voltage 

a 

across the device. Since the differential change in gate 
charge is equal to the insulator capacitance rimes the 
differential change in voltage across the insulator, the 
capacitance can be written as 



C = 



C. 

1 




( 17 ) 



where Vg is the voltage across the silicon surface region. 
By rearranging the terms and integrating, this equation 
becomes 



V 



S 





dV^ + K 

d 



= Ug + K, where is the integral. 



( 18 ) 



If the high frequency normalized capacitance is considered, 
so as to neglect the contribution from interface 
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capacitance C„„> Uq becomes 

o o ^ 



“s ' 



'c.’ 

1 



dV 



+ 1 



(19) 



The quantity (C./C„)'^ can be obtained from the experimen- 

X o 

tal high-frequency data from the relation 



'C. 
1 



'\ 2 






c. 

1 



- 1 



HP 



( 20 ) 



The (C^/Cg)^ - Ug plot is obtained by evaluating Eqs. 19 

and 20 and is the first plot drawn by the computer. 

A least-squared-error fit is applied to the linear 

l6 

portion of this plot, if it exists. From the slope of 
this line the surface impurity concentration can be found by 

r. 2 



Ns = - 2 



£ qA m 



( 21 ) 



where m is the slope. In the above equations, the surface 
was assumed to be in depletion" so the surface potential 



The linearity of the (C^/Cg)^ versus Ug plot is an 
indication of two important properties of the MIS. If the 
curve is not linear or is linear in two regions it may in- 
dicate an inordinate amount of interface states or a gross 
nonuniformity in the insulator or semiconductor. 

17 

Depletion region markers showing the depletion region 
used are plotted on the C-V high-frequency curve output and 
determined by the criterion (.050^ + .95Cmin) -Cl (.25Ci+ 

.75Cuij_n) where Cmin is the minimum capacitance value. 
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could be considered zero and the Inverse of the space charge 
capacitance versus V plot has a constant slope. 

The impurity concentration in the bulk can be cal- 
culated from the following equation which is valid when the 
surface capacitance is at its minimum value. 




n^exp 



2kT J 



■ 2 





( 22 ) 



Cgmin obtained from the experimental high frequency 

capacitance curve and is related to Vg by 



fc 



kT 



0 



Smin 



A 



exp [q 



Vr 



/2kT] 



[2(q|V^ /kT - 2)] 



(23) 



With calculated and using the relationship 



'*>b 





(24) 



where is the potential in the bulk, the constant, K in 
Eq. 18, can be determined and Vg obtained as a function of 
V„ . This is the third plot output of the program. The 
program also provides tabular data of Ng, and 4>^. 

2 . Interface/State Density 

The Interface state density can be found from the 
following equation 



N = — r 
SS qA ^SS 

where Ngg is given in units of 
found from 



(25) 



[volt -meters ^ ^ and Cgg is 
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m 









c. 

1 



( 26 ) 



'SS 



fc. 1 




fC. I 


1 




1 


c 


*“ ± 
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w J 


LF 


J 



- 1 



HF 



The computer solution to Eq. 25 Is the fourth graph output 
of the program with the horizontal axis being defined as 
-(<j>b is valid only in u’r interval -- 24 >b<Vg <0 when 

the semiconductor is in depletion. Tine progz^am also outputs, 
in tabular form, three values of Xgg corresponding to Vg 
values of approximately - 24*^3 ”9^ 

3 • Effective Interface Charge 

The effective interface charge is a form of insula- 
tor charge and is given by the equation 



Q 



eff 




(27) 



where Q can be determined from experimental data as follows 
s 



Qg 



fc 



- 1 



HF 



- N, 



exp 



kT 



+ 2 



N. 

1 



N^ 



qVc 



sinh — 



kT 



( 28 ) 



Equation 27 is the fifth plot output from the program. 
Tabular data for 3 values of and are provided for 

Vg values of ” 24)^3 -4)^ and 0. 
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B. COMPUTER PROGRAM BASED ON HIGH-FREQUENCY C-V DATA ONLY 



There are several papers on the evaluation of MIS 
structures from high-frequency curves Cl5, l8, 30] which 
provides the following information: 1) conductivity type, 

2) impurity concentration, 3) minority carrier lifetime and 
4) interface state density. Quantities 1, 2, and ^ are 
obtainable directly from C-V curves. Quantity 3 is ob- 
tained from the C-t curve when a MIS is biased rapidly from 
accumulation to inversion. Other quantities obtainable 
from C-V curves at different temperatures and conditions 
are charge transport in the Insulator, trapping effects, 

polarization effects, and space charge build-up due to ionlz- 
1 8 

Ing radiation. Except for the last effect these quantities 
will not be discussed in this thesis. 

1 . Conductivity Type 

The determination of semiconductor conductivity 
type is a simple matter from a high-frequency plot. One 
simply observes the slope of the C-V curve as the dc bias 
is swept from accumulation to inversion (or vice versa). 

If the slope is negative, the semiconductor is p-type and 
if the slope is positive it is n-type. This method for 
determining conductivity type is faster than the Hall 



1 8 

The ionizing radiation effect is discussed in Chapter 
IX of this thesis. In that chapter C-V curves are evaluated 
before and after electron Irradiations of a standard CN I 38 
N-channel depletion MOSFET. 
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19 

method, however, a drawback is the required insulator 
film deposition. A film deposition is not, of course, 
required in the Hall method. 

2 . Interface State Density 

When the deviation from the ideal C-V curve is the 
result of surface states, a determination of these states 
can be found by the amount of shifting that takes place on 
the voltage axis of the C-V plot. A computer program that 
first calculates the ideal C-V curve and then compares it 
with input experimental C-V curves is discussed in t?iis 
section [ 30 ]* 

The following input parameters are used: the accumu- 

lation capacitance in picofarads, the impurity concentration 
in cm , the temperature in degrees Kelvin, the composition 
X or Y (for lead-tin semiconductors), the dielectric con- 
stants for both the semiconductor and insulator, the insula- 
tor thickness in Angstroms and the experimental C-V in 

20 

either normalized or absolute quantities. 



19 

The Hall method has been the standard procedure for 
determining conductivity type. It entails applying a mag- 
netic field across the semiconductor in which a current is 
flowing and measuring the Hall voltage which is set up 
because of carrier drift in the magnetic field. This method 
requires samples in the form of Hall bars and reduced tem- 
perature. A by-product, perhaps the major product, of the 
Hall measurement is the impurity concentration. 

20 

Capacitance values in the high-frequency Inversion 
region must be input as negative quantities. 
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The following quantities are the outputs of the 
program: the Fermi level above the valence band, the nor- 

malized high and low-frequency inversion m.lnlmum capacl- 
21 

tance, values of capacitance, values of the surface 

potential in volts, the energy above the valence band to 

which the interface states are filled in ev., and the 

_2 

density of interface states in cm 

This program was used in this thesis to evaluate 
the cases of n-type Ge-MIS and the lead-tin alloy MIS 
samples whose low-frequency C-V could not be measured. 



C. COMPUTER PROGRAM TO ANALYZE HIGH-FREQUENCY INVERSION 
MINIMUM CAPACITANCE 

In the high-frequency inversion C-V case a minimum 
capacitance exists and is defined by the equation 



C . 



c. 

1 



mlnHF C. 



^Slnv 



+ 1 



where C_ . 

Sinv 



""s S 

^inv 



^inv 



f 4Kg kT In (N/n. )'» 
q^N 



( 29 ) 



21 

These values are presented to allow a comparison of 
actual versus theoretical minimum capacitance values. 
Their agreement depends, to a great extent, on the value 
of impurity concentration input. 
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In order to determine the combinations of K. , K_ , N 
and which will produce the low experimental values of 
CminHp lead-tin samples, a computer program was 

developed to calculate ^.s a function of impurity 

concentration (N), insulator dielectric constant (K^. ) , 
semiconductor dielectric constant (Kg) and insulator thick- 
ness (X^). The result of this analysis appears in Chapter 



VII. 
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VI. PREPARATION OP MIS SAMPLES OP IV -VI SEMICONDUCTORS 



The device fabrication procedure is extremely critical. 

In the well developed silicon industry, it is relatively 

easy to grow or deposit an insulating layer over a silicon 

wafer and then vacuum deposit the metal layer. However, 

even with silicon, it is still difficult to produce repeat- 

able MOS devices. When applied to IV-VI semiconducting 

materials, this difficulty is more severe since the MIS 

properties of these elements have never been studied. 

Methods for applying a good insulating film are not known. 

A good candidate is Rutile (TiOj) whose dielectric constant 

22 

is commensurate with the IV-VI semiconductors. Presently 
Dr. N. Bottka at the Naval Weapons Center (NWC), China Lake 
is attempting to develop a deposition procedure for TiO^. 
This section of the thesis contains the fabrication pro- 
cedure for construction of the MIS samples. 

A. Pb^_^Snj^Te AND Pb^_YSnYSe THIN PILM SEMICONDUCTORS 

The Pb^ ^Sn^Te(Se) studied in this research was in the 
form of thin films deposited on CaPj, BaPj, and KCL substra- 
tes using the one-boat deposition method and the Knudson 



22 

The lead-tin sem.iccnductors have a dielectric constant 
in the range of from 100 to 400. TiO 2 has a reported dle- 
elctric constant up to 100. Dielectric constants of 40 have 
been reported at Western Electric and a stable 80 at the 
University of New Mexico [12]. 
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cell deposition method. Metallurgical measurements were 
then made on the samples in order to determine their crystal 
structure, orientation, and thickness. Conventional elec- 
trical measurements based on the Hall effect were then made 
from 300 ° K to 90° K to determine the carrier type, carrier 
concentration and Hall mobility. 

The sample preparation, metallurgical evaluation and 
electrical measurements were not performed by the author 
and will not be described in this thesis. A complete des- 
cription of these processes can be found in Ref. 11. Both 
single crystal and polycrystalline thin films were used. 

The mobilities of the thin films were in the range of from 
80 to 29,000 cm^-v ^-sec ^ and the concentrations (p and n) 
ranged from 10*® to 10*® cm Figure 6-1 gives the shape 
of the samples used; one a Hall bar for Hall measurements, 
the other an optical sample for optical measurements. 

B. INSULATORS 

Insulator deposition is the most critical step in pro- 
ducing useful MIS structures. A poor MIS causes C-V anoma- 
lities such as drift with time, non-repeat able 
characteristics, and low-voltage breakdown along with the 
possibility of having no capacitance change at all. Causes 
of these anomalltles are pinholes through the insulator, 
non-uniform insulator thickness, mobile ions in insulator 
and Interface states. 

Si-MOS fabrication techniques are well established. The 
major Insulator material used is SiOj and is produced by the 
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INSULATOR 





Figure 6-1. MIS Thin-Film Sample 
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oxidation of the surface of the semiconductor. For less 
well developed semiconductors such as the III-V and IV-VI 
semiconductors, other techniques had to be developed. 
Extensive MIS fabrication procedural work has been conducted 
at the University of New Mexico for GaAsP samples. Various 
insulator deposition m.ethods such as RF sputtering of SIO 2 , 
Al^Oj, SijN^ and TiO^ and wet oxidation were tried and dis- 
carded after they were found to produce unstable insulators. 
The technique that proved to be most successful was the 
passing of atmospheric-pressure-dry oxygen over the sample 
at temperatures of 710° C or lower and then annealing the 
samples in nitrogen for more than 40 minutes. The cleanli- 
ness of the sample is critical. Mechanical lapping and 
polishing accompanied by chemical cleaning and polishing 
of the semiconductor surface before insulator deposition is 
an involved and complicated process [19]- 

The deposition scheme used for the Pb^_^Sn^Te (or 
Pbi_YSnYSe) samples consisted of first cleaning the samples 
with trichloroethyne , then ethanol and finally a rinsing 
in heated deionized water. The samples were placed in small 
flasks of trichloroethyne and hand agitated for 3 minutes. 
The trichloroethyne was replaced in the flask by ethanol and 
agitated again for 3 minutes. The final cleaning step con- 
sisted of replacing the ethanol by heated deionized water 
and agitating before the samples were removed and dried by 
passing dry N 2 over them. Metal masks were then put over 



78 



the samples and loaded into the vacuum system as soon as 
possible to avoid surface contamination. 

The deposition equipment consisted of a Norton vacuum 
station with a Temescal 15 KV. E-gun and Sloan deposition 
thickness monitor model DTM-3. The system vxas evacuated 

-7 

overnight, to approximately 10 Torr. The E-gun current 
was next increased slowly to out gas the evaporant and pre- 
vent sputtering. 

The final step before the actual evaporation was to fo- 
cus the E-beam spot size on the evaporant to ensure proper 
heating. With the shutter in place to prevent sporadic 
evaporation until the evaporant was heated properly, the 
E-gun was brought up to operating pov/er v/hlch was approxi- 
mately 20^ of the maximum 15 KW full power. When the evap- 
orant was observed to be properly melted and the Sloan 
model DTM-3 deposition thickness monitor was set to its 
zero level, the shutter was removed and the deposition 
began. The Sloan meter frequency was observed, a stop watch 

was started and evaporation continued for a predetei’mined 

23 

frequency or tim.e. The sample v;as then allowed to cool to 



23 

•^The principle of operation of the DTM-3 is to mix the 
frequency of a stable local oscillator in the meter with the 
frequency of a crystal oscillator in the deposition jar, 
whose frequency is dependent on the amount of mass on it. 

The beat frequency of these two oscillators determines the 
thickness of the film deposited by the following formula: 

where : 

_ 2Af T = thickness in Angstrons 

p Af = beat frequency 

p = density of deposited material 
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room temperature which took typically 1 hour and the films 
were checked for physical characteristics. Table I 
provides a list of the samples prepared and the materials 
used for insulator. Adhesion properties on all samples were 
good except for ZnS. The insulator surface appeared to con- 
tain a fair but not inordinate amount of obvious pinholes. 



Semiconductor 
Pbi__^Sn^Te (Se) 

rt 

M 

ft 

tt 

rr 

Ge (p) 

Ge (n) 

Si 

Table I. List of 



Insulator 

AI2O3 

AI2O3 

AI 2 O 3 

Si02 

ZnS 

Ti02 

AI 2 O 3 

AI2O3 

SIO 2 

Samples Prepared 



Thickness (A°) 

100 

300 

980 

450 

1000^^ 

500^5 

250 

180 

1000^6 

:• Experimentation 



ZnS was deposited by thermal heating in the Physics 
Laboratory at the Naval Postgraduate School with the help of 
Mr. Donald Spiel. 

25 

Ti 02 has not yet been deposited; however, the deposi- 
tion technique is being developed by Dr. Botka at NWC, 

China Lake. 



The SiOj over silicon was grown at Stanford University. 
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C . METAL DEPOSITION 



Metal deposition was in the form of Ni pads thermally 
deposited through various shaped masks. The various shaped 
masks afforded the opportunity to observe the effect of 
different gate areas on experimental curves. Deposition was 
performed similar to the above procedure. Depending on the 
thickness of evaporant required, in this case two small one- 
quarter inch pieces of 50 mil wire, the evaporant v;as put in 
a boat and thermally heated first to form the two pieces 
into a single bead and to out gas and then finally to evap- 
orate when the proper vacuum of lO”^ was reached. Thickness 
was monitored by measuring the evaporant resistance across 
a 2 by 5 cm area. V/hen this resistance reached 1 Kfi,the 
evaporation process was terminated. In all cases the 
evaporated metal adhered properly. In several cases the 
metal was barely visible because of its thinness. 

Three other type metal gates were used. On the Si 

sample 30 mil. and 20 mil diameter aluminum gates were 

deposited. On other samples no metal deposition was made 

27 

and a small drop of Hg was put on the tip of a probe to 
act as a metal gate. The Hg probe was also used on several 
samples which had other Ni metalization which because of the 
thinness of the sample, allowed the point probe to punch 
through the metal. The third type of metallization was in 



27 

Here there was a problem defining the area accurately 
since the area depended on the force applied to the 
mechanical probe. 
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the form of Au gates of 30 mil. diameter and on several 
Hall samples in the form of strips covering the Hall bars. 
The gold deposition did not adhere well. The Hg probe has 
given more repeatable results. 
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VII. EXPERIMENTAL RESULTS OF 



IV-VI ALLOY SEMICONDUCTOR MIS 

The experimental results of IV-VI alloy semiconductor 
MIS are presented in this chapter. These results consist 
primarily of high-frequency C-V curves. The low-frequency 
and C-t response were considered, attempted but not pre- 
sented, except in one case for each, since successful plots 
were not obtained. Several G-V plots were also attempted. 

A. MIS SAMPLES STUDIED 

MIS samples of 6 different IV-VI semiconductors were 
experimentally investigated. Two different types of insula- 
tors of three thicknesses were used. These combinations 
were: PbTe with Al20^ and Si02, Pb is"^® with Al202» 

Pb QgSn 20'^^ with AI2O2 Si02, Pb ygSn 24'^® with Al202> 

PbSe and Pb ^ySn with AI2O2, and Pb ggSn with 

Si02* They are also listed in Tables II and III. The major 
portion of the measurements were made at room temperature. 
Several samples were tested at liquid nitrogen temperature. 

Since most of the samples showed leakage through the 
insulator, the measurement method using the lock-in-ampli- 
fler was not useable and, consequently, only a few low- 
frequency C-V measurements have been obtained. G-V plots 
were also attempted on several samples but only one success- 
ful curve was obtained. 
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B . HIGH-FREQUENCY C-V RESULTS 



This section contains the experimental high-frequency 
C-V measurements at room temperature and at liquid nitrogen 
temperature and some low-frequency data. A theoretical C-V 
curve is first presented and followed by the experimental 
curves . 

1. Pb^_^Sn^Te MIS C-V curves at Room Temperature 

The theoretical C-V curves for the PbTe and Pb 
Pb^ samples at 300°K is presented in Figs. 7-1 

and 7-2. The material parameters selected are the following: 
impurity concentration assumed was 1X10*® cm semicond- 

uctor dielectric constant assumed was 400. The high frequen- 
cy C-V experimental data are presented as follows. 



MIS 

Semiconductor 


MIS Insulator 
and Thickness (A°) 


Figure Number 


PbTe 


AI 2 O 3 ( 100 ) 


7-3 




SIO 2 (450) 


7-4 


^^. 82 ^*^. 18 ^® 


AI 2 O 3 ( 300 ) 


7-5 




AI 2 O, ( 100 ) 


7-6 




AI2O3 (300) 


7-7 




SiO 2 (^ 50 ) 


7-8 


^^. 76 ^^. 24^® 


AI 2 O 3 ( 100 ) 


7-9 




AI 2 O 3 (300) 


7-10 



Table II. Pb^-x^nxTe Semiconductor MIS Experimentally 
Investigated 
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Figure 7-1. Theoretical C-V Curve for PbTe MIS (Ideal), N = 1X10*® cm 
T = 300°K, Insulator: Al^O, (100A°) 
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Figure 7-2. Theoretical C^^V Curve for Ph ^ 00^^ . 20*^® (Ideal) 

N = 1X10^® cm~^, T = 300°K, Insulator: S102(^50A°) 
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Figure 7-'^. Experimental C-V Curve for PbTe MIS, T = 300°K, Insulator: 
SIO, (450A°) 
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Figure 7-5* Experimental C-V Curve for Pb gpSn i sTe MIS, T ?00^K 
Insulator: AIjO, (300A°) 
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Figure 7-6. Experimental C-V Curve for Pb 20'^® MIS, T = 300°K 

Insulator: A1,0, (100A°) 
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Figure 7-7. Experimental C-V Curve for Pb ^ 20*^® MIS, T = 300°K 

Insulator: Al.O, (300A°) 
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Figure 7 - 8 . Experimental C-V Curve for Pb_ 80 ^^. 20 Te MIS, T = 300°K 
Insulator: Si02 (^50A°) 
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Figure 7 - 9 . Experimental C-V Curve for Pb . ygSn . 2 i|Te MIS, T = 300°K 
Insulator: A1.0_ (100A°) 
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Figure 7-10. Experimental C-V Curve for Pb 7 gSn 24"^® MIS, T = 300°K 
Insulator: A1,0, (300A®) 



There are several differences betv;een theoretical 



and experimental curves. First j the measured high-frequency 
minimum capacitance is much lower than the theoretical 
value. Possible explanations for this capacitance differ- 
ence are the following; nonuniformities in the insulator 
or semiconductor, the semiconductor dielectric constant may 
not be as high as 400, error of the impedance measurement in 
the determination of capacitance brought out by the leaky 
insulator. The latest factor was discussed in detail in 
Chapter V and is, in most probability, the major cause for 
the differences. Although this difficulty with leaky insula- 
tors exists and the curves are distorted, depletion and in- 
version layers v;hich were controlled by the gate voltage 
seem evident from these experimental results. 

Second, drift was common in almost all the samples. 
Drift is not Included in these figures except in Fig. 7-6 
where it is presented for gate 3* Third, hysteresis ef feres 
sometimes were found. Their effects were generally not 
large and are not presented. 

Fourth, there is often a small hump in the C-V curve 
near zero as shown in Figs. 7-3 j 7-4, 7-8, and 7-10. 

Fifth, the voltage bias ranges from the accumulation 
to Inversion is much larger than the theoretical results. 

i 

Sixth, low-frequency and high-frequency deep deple- 
tion C-V curves can not be obtained. 

Figure 7-9 shows a G-V curve in addition to the high- 
frequency C-V of a sample. G-V measuremnnts were not 
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attempted in general, since G-V curves do not provide as 
much information about the sample MIS as C-V curves do. 

Future success v;ith G-C curves could prove to be of value if 
lovj-frequency curves are not possible. Errors associated 
with phase-setting on the LIA G-V curves for leaky insulator 
MIS is not as critical as it is in C-V curves. This was 
discussed in Chapter V. 

2 . Pb _ ^ Sn^Te MIS C-V Curves at Liquid Nitrogen 

Temperature 

Considerable difficulty was experienced in obtaining 
low-temperature C-V curves. First, Hg freezes at liquid 
nitrogen temperature and the measurements become irratic in 
many cases. Second, the samples were immersed in the liquid 
nitrogen in the cold thermos and the bubbling of the liquid 
nitrogen near the sample caused jitter in the measurement. 

Figure 7-11 shows the temperature variation of the 
high-frequency C-V of a PbTe ( A1 2 O 3 )MIS . The lowest tempera- 
ture is labeled T1 in the figure. The accumulation capaci- 
tance was found to decrease with temperature. It contradl<. ts 
the theory which asserts that accumulation capacitance re- 
mains constant for all temperatures. Instead, the inversion 
capacitance should decrease with temperatures since the 
minority carrier lifetime is reduced. This abnormality of 
temperature variation was not characteristic of all the 
samples. Another sample. Fig. 7-12, for example, shows two 
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Figure 7-11. Experimental C-V Curve for PbTe MIS at Various Low Temperatures 
Insulator: Al-O, (300A°) 
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Figure 7-12. Experimental C-V Curve for PbTe MIS, T = 80°K, T = 300°K 
for Gate 1, Insulator: Al.O, (100A°) 



In this 



low-temperature curves for different gate areas, 
case the trend was the opposite and the inversion capaci- 
tance went down with temperature. This fact is not clear in 
Fig. 7-12 since the capacitance is normalized to the insula- 
tor capacitance. The absolute capacitance values in this 
case were 6 picofarads in accumulation and 2 picofarads in 
inversion. The accumulation capacitance represents a dras- 
tic reduction of the typical measured accumulation capaci- 
tance at room temperature. The major difference present in 
these curves over the room temperature case was the percen- 
tage change between accumulation and inversion capacitance. 
In the majority of the room temperature cases, the percen- 
tage change was ^ 0 %. The reason for this difference is un- 
known. In general, low-temperature curves have been 
difficult to obtain. Since the lead-tin semiconductor MIS 
will be operated at cooled temperature, low-temperature 
measurement must be determined in the future. 

3 . Pb ^ y SUySe MIS C-V Curves at Room Temperature 
The theoretical C-V curves at 300°K for a 
PbQ 1^^ (Si02~^50A°) is presented in Fig. 7-l3- The 

carrier concentration was 1.0X10^® cm The high-frequency 
experimental C-V curves are presented as follows: 



? 8 

°Gate 1 in Fig. 7-12 represents the room temperature, 
case for the sample. 
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Figure 7-13. Theoretical C-V Curve for PbQ oSno iSe MIS (Ideal), 
N = 1X10^® cm-3, T = 300°K, Insulator: SIO (450A°) 
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MIS 

Semiconductor 


MIS Insulator 
and Thickness (A°) 


Figure Number 


PbSe 


AI 2 O 3 


( 300 ) 


1 — 1 

1 




AI2O3 


( 300 ) 


7-15 


^'^.97®". 03®® 


A 1,0 3 


( 300 ) 


7-16 




Si 02 


(450) 


7-17 



Table III. Pbi-ySnySe Semiconductor MIS Experimentally 
Investigated 

The difference between the theoretical and experimental 
curves is obvious from Figs. 7-13 and 7-1^. These differ- 
ences are the same as those discussed for PbSnTe; however, 
they are even larger because of the small capacitance varia- 
tion in the theoretical curve. All experimental C-V curves 
contain a large capacitance change from accumulation to in- 
version. The only low-frequency C-V curve obtainable ap- 
pears in Fig. 7-15 for a PbSe (AI 2 O 3 -300A° ) sample. The 
low-frequency was obtained at lOKHz , but it is not the tru' 
low-frequency curve since the Inversion and accumulation 
capacitances are not equal. Lower frequency curves were 
tried but not obtainable. 

4 . Pb ,^ Sn, ^, Se C-t Response 

The capacitance time response was attempted. But, 
because of the short minority carrier lifetime of the samples 
at room temperature, no results were measurable except for 
a few cases. Figure 7-18 gives the C-t response for 
Pb ggSn Q^gSe pulsed from -15. V to O.V at room temperature. 
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Figure 7-1^* Experimental C-V Curves for PbSe #1, T = 300°K, Insulator 
A1,0, (300A'’) 
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Figure 7-15. Experimental C-V Curves for PbSe MIS #2 Showing Low- 

Frequency Response, T = 300°K, Insulator: A1 0 (300A°) 



( 




10^1 



Figure 7-16. Experimental C-V Curve for Pb qySn oS^e MIS, T = 300°K 
Insulator: A1,0, (300A°) 
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Figure 7-17* Experimental C-V Curves for Pb goSn loSe MIS, T = 300°K 
Insulator: SIO, (450A°) 
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Figure 7-18. C-t Response for Pb i^Se MIS, T = 300°K, Insulator: 

SiO, (i^50A°) *5° ■ ^ 



i 



It is clear from the figure that the time for the MIS to 
relax from the deep depletion to the quasi-equilibrium in- 
version is about 0.1 sec. This storage time of 0.1 sec. 
is much longer than the expected time which is on the order 
of a millisecond or less. The validity of this measurement 
is questioned; however, the same result was obtained at 
various different gate locations on the sample. Storage 
times of this magnitude would allow normal integration and 
readout times for an IRCCD application. 

C. ANALYSIS OP MEASURED C-V DATA 

Because low-frequency results were not satisfactory for 
the lead-tin semiconductor MIS samples, the computer pro- 
gram developed for use with both high and low-frequency C-V 
curves could not be used except on one sample where a low- 
frequency C-V curve was obtained. The program based on high- 
frequency C-V data only, was also tried. It too container 
errors because the minimum experimental values of capaci- 
tance were much lower than theoretical values. The thlrc 
part of this section contains computer calculations of 
minimum high-frequency capacitance versus Impurity concen- 
tration with insulator type and thickness as parameters. 

The purpose of this problem is to determine what combinations 
of MIS parameters of the semiconductor and insulator would 
lead to a theoretical C-V in better agreement with the measured C-V results. 

1 . Low and High-Prequency Data Analysis 

The data analysis com.puter program discussed in 
Chapter V was applied to the C-V curves of the PbSe sample. 
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Fig. 7-15 • Because the curves are distorted from normal 
C-V variation, it was expected that abnormal behavior 
would result. 

No useful output curves were obtained. The tabular 
data output appears in Table IV. The Impurity concentration 



Quantity 


Value 


ENS 


-1.62X10^® cm~^ 


ENB 


-1.37x10^^ cm“^ 


*^b 


-0.l89v 


N (at - -0.083v) 

ss 


1.24X10^2 cm"^ 


N (at ip = O.Ov) 

Q 0 • 


9.35x10“ cm"^ 



Table IV. Summary of Data Analysis Output for PbSe 

(AI 2 O 3-100A° ) MIS from Experimental High and Low- 
Frequency C-V. 

does not agree with expected values. For exam.ple, the 

impurity concentration in the bulk (ENB) is less than the 

intrinsic concentration (about 10^^ cm~^). The N value i . 

s s 

about the value expected. 

2 . High-Frequency Only Data Analysis 

The high-frequency only data analysis described in 
Chapter V was applied to the high-frequency curves of four 
samples: PbTe (S102-450A° ) , Pb qqSu (Si02-450A° ) , 

Pb ^gSn 2 i|Te (Al^O j-lOOA® ) , and PbSe (AlzOa-SOOA® ) . The 
results of this analysis appear in Table V. Table V shows 
that the normalized theoretical minimum high-frequency 
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Sample Bias (V) Normalized Surface Potential Interface States 

^MINHF (ev) (cm~^) 
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Semiconductor MIS Samples. 



capacitance is much larger than that obtained experimentally 

therefore, errors are expected from this data analysis. For 

the Pb ygSn 24 '^® sample the theoretical minimum capacitances 

is lower by a considerable amount and more accurate results 

were obtained. Note that the interface state density in 

_2 

all the samples was about 10^^ cm . This density is 2 to 
3 orders of magnitude higher than that obtained from good 
quality Si-MOS. 

3 . Theoretical Minimum High-Frequency Capacitance 

Analysis 

Because of the large difference between theoretical 
and experimental high-frequency minimum capacitances, the 
program, discussed in Chapter V, was developed to determine 
what combinations of insulator thickness. Insulator type, 
semiconductor dielectric constant and impurity concentra- 
tion would lead to better agreement of these capacitance 
values. Each set of data pertains to a certain insulator 
type and thickness and the minimum capacitance versus im- 
purity concentration is plotted with the semlbonductor 
dielectric constant as a parameter. The results of this 
analysis appear in Fig. 7-19 for AljOj (100A°) and in Fig. 
7-20 for Ti 02 (100A°). These two figures represent the 
closest agreement with the theoretical minimum capacitance. 
Table VI gives the absolute capacitance minima. Insulator 
TiO^ was considered because its high dielectric constant 
makes it compatable with the lead-tin semiconductor. Table 
VI indicates that TiOj at 1,000A° gives about the same 
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Figure 7-19. Plots of Cjyijj^j^p Versus Impurity Concentration for IV-VI 
Semiconductor - A1 0 (100A°) MIS 
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Figure 7-20. Plots of C|vjj«,up Versus Impurity Concentration for IV-VI 
Semiconductor TiOj (100A°) 
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Impurity 



Insulator Semiconductor Minimum Value Concentration 



Type 


Dielectric 

Constant 


of Normalized 
‘^MINHP 


(cm~3) at 
Minimum 


AI 2 O 3 (100A°) 


100 


0 .il 6 l 


4X10*® 


AI 2 O 3 (300A°) 


100 


0.720 


4X10*® 


Si02 (450A°) 


100 


0.897 


4X10*® 


TIO 2 (100A°) 


100 


0.086 


4X10*® 


Ti02 (500A°) 


100 


0.320 


4X10*® 


Ti02 (1,000A° 


) 100 


0.488 


4X10*® 



Table VI. Summary of Absolute Capacitance Minima for IV-VI 
Semiconductor MIS 



value of minimum capacitance as does 100A° of AI 2 O 3 . The 
best fit with experimental minimum capacitance is 100A° 
TiOj . 
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VIII. EXPERIMENTAL RESULTS OF GERMANIUM MIS 



Experimental analysis for two different types of Ge-MIS 
is described in this chapter. To date, almost no research 
has been done on Ge-MIS. Si-MOS is v;ell established and one 
doesn't gain much of an advantage from Ge with its energy gap 
of 0.67 ev. , except that its spectral response moves farther 
toward the infrared. The study of Ge-MIS by the author is 
a joint effort with Dr. N. Bottka from NWC , China Lake who 
is doing electro-reflectance research for both basic band 
structure study and applied tunable IR filter development. 

A. DEVICE STRUCTURE 

The devices studied were a p-type Ge-MIS with a resisti- 
vity of 0.05fi-cm and an n-type Ge-MIS with a resistivity of 
I 29 

40Q-cm. The samples were rectangular in shape of dimen- 
sions 7 by 17 mm. The p-type sample had 250A° of Al^O^ anc 
the n-type had a l80A° of Al^O^ as insulator, both depositf i 
by E-gun evaporation techniques at NWC, China Lake. The 3 
by 9 mm Ni metal gates were deposited by resistive heating. 
Copper metal electrodes were epoxyed to both ends of the 
metal gate and shellac was put over the surface of the gate. 
These electrodes did not provide satisfactory results for 



29 

These values of resistivity correspond to impurity 
concentrations of 1.3X10^^ cm~3 and 5.0X10^^ cm~3 respectively. 



MIS studies. On one sample an electrode did not make good 
contact with the metal gate and no experimental results 
were possible from it. The metal electrodes on both samples 
broke and the Hg probe method was finally used. 

B. THEORETICAL PROCEDURE 

The theoretical procedure for analyzing the Ge-MIS con- 
sisted of obtaining first the ideal low-frequency inversion, 
high-frequency inversion and high-frequency deep depletion 
ideal C-V curves. These curves were generated by the program 
developed by Doshier [8]. The calculated C-V curves appear 
in Pig. 8-1 for n-type Ge-MIS and in Fig. 8-2 for p-type 
Ge-MIS. The tabular output of the computer program for 
p-type Ge-MIS appears on page l60. 

The p-type Ge-MIS C-V plot is typical of the curves des- 
cribed in Chapter III of this thesis. Q storage is obtained 
from the difference between the high-frequency inversion and 
deep depletion capacitance at a certain gate voltage. A 
typical value at = 7* 53V, as seen on page l60, is 3*04} ^0^^ 

u 

-2 

coul-cm , which is comparable to Q storage for Si-MOS. 

Other quantifies shown on page l60 are the depletion and 
deep depletion depths at various gate voltages, the flat- 
band capacitance (0.774) and the surface potential at the 
onset of inversion (0.449V). 

The n-type Ge-MIS C-V plot is not typical of the curves 
described in Chapter III. The primary reason for the 
difference is due to the low impurity concentration of the 
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Figure 8-1. Theoretical C-V Curves for n-Type Ge-MIS (Ideal), Insulator: 
A1,0, - 180A°, N = 5X10 cra-3, T = 300°K 
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Figure 8-2. Theoretical C-V Curves for p-Type Ge-MIS (Ideal), Insulator: 
A1,0, - 250A°, N = 1.3X10^’ cm-3 T = 300°K 



sample, was 2.2X10“* cm ^ whereas the impurity concentra- 

13 

tion vms 5X10 cm . This sample would not be useful as 
a CCD. Because its high frequency depletion and inversion 
capacitances are nearly equal, the storage capacity is very 
small . 

C . EXPERIMENTAL RESULTS 

The experimental procedure consisted of three parts: 

1) determination of high and low-frequency C-V curves, 

2) C-t measurement and 3) the effect of shunt resistance on 
the experimental curves. This experimentation was conducted 
as described in Chapter IV on experimental methods. 

1 . High and Low-Frequency C-V 
a. P-type Ge-MIS 

The experimental p-type Ge-MIS C-V curves agreed 
generally with the theoretical curves with three differences 
as shown in Fig. 8-3. First, a hysteresis effect was found. 
Second, the measured minimum high-frequency capacitance is 
larger than the theoretical calculation. Third, the fre- 
quency used in measuring the low-frequency C-V is not low 
enough because the normalized low-frequency inversion 
capacitance is less than one. 

The C-V data of this p-type Ge-MIS sample is 
not uniform as shown in Fig. 8-4 for three different gate 
locations. The differences in normalized C are evident. 

In section D, they will be analyzed separately. Their low- 
frequency curves are also different but not shown in this 
figure . 
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Figure 8-3. Experimental C-V Curves for p-Type Ge-MIS Showing the Hysteresis 
Effect, Insulator: A1,0, - 250A®, N = 1 . 3X10 cm"3 , T = 300°K 
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Figure 8-4. Experimental C-V Curves for p-Type Ge-MIS Showing Nonuniformi 
at Different Gate Locations, Insulator: AljOg - 250A°, N = 
5X10“ cm-3, T = 300°K 
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b. N-Type Ge-MIS 



Figure 8-5 gives the C-V curves of the n-type 

Ge-MIS. They differ from the theoretical calculation in 

two ways. First, the measured inversion capacitances are 

higher than calculated values. This can be attributed to 

high contact resistance . Second, the gate voltage range 

over which the accumulation to inversion capacitance change 

is much larger. Experimentally this range v;as between -8.0 

and l6.0 volts whereas theoretically this, range was between 

0.34 and -O .56 volts. A possible explanation of this dif- 

31 

ference is an over abundance of interface states. The 
experimental curves are atypical of standard Si-MOS curves. 
Although hysteresis also exists as shown in Fig. 8-5, it 
was not as noticeable as the hysteresis for the p-type 
Ge-MIS. 

Since this n-type Ge is almost intrinsic, both 
the dark current and the generation of minority carriers 
are small. Consequently, low-frequency C-V variation was 
not observed. ‘ 

2 . C-t Measurement 

Figure 8-6 presents the experimental C-t response 
of the p-type sample, along with the high-frequency C-V. 

The voltage pulse was from -8. v to +8.v. The small sharp 



^^"Finger print” Fig. 3-6b shows this result. 

^^The high-frequency analysis in section D.2 shows about 
5 times the Interface state density over the p-type Ge-MIS. 
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Figure 8-5. Experimental C-V Curves for n-Type Ge-MIS, Red: Reduced Scale 
Shovring Hysteresis, Exo: Expanded Scale, Insulator: AI2O3 - 
180A°, N = 5 x 10 *® cm-3,* T = 300°K 
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Figure 8-6. Experimental High-Frequency C-V and C-t Response for Ge-MIS 
(p-Type), a) C-V Response, b) C-t Response 
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spike at t = 0 sec is probably the result of the filling 
of fast interface states. The relaxation time (or storage 
time) is very long as noted by the gradual rise of the 
capacitance. In this case, it was approximately 60 sec. 

This is the time for the capacitance to attain the high- 
frequency inversion capacitance from the deep depletion 
capacitance. No C-t response could be measured on the 
n-type Ge-MIS. 

D. ANALYSIS OP Ge-MIS DATA 

Data analysis consisted of evaluating the low and high- 
frequency C-V data using the computer program described in 
Chapter V, section A. For the p-type sample three sets of 
data, from the same sample, were analyzed. Using the high- 
frequency C-V, both the p-type and n-type Ge-MIS were 
analyzed by the program described in Chapter V, section B. 

1 . P-type Data Analysis 

The results for one gate on the p-type Ge-MIS appear 
in Figs. 8-7 through 8-11. Figure 8-7 is the absolute low 
and high-frequency C-V curves. Figure 8-8 gives (C,/C„)^ 

J. O 

as a function Ug together with the straight line fit to 
the curve. Figure 8-9 is a plot of V„ versus V . Figure 

o 3. 

8-10 shows the relationship of the interface state density, 
Ngg, and the surface potential. Figure 8-11 provides a 
plot of versus Vg. 

An examination of these curves reveals that the 
flat-band voltage, where Vg equals zero in the V^-V^^ plot. 
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Is -1.2 volts and that the range of values of Vg extends 

from -0.35 volts to +0.45 volts. The total range of Vg 

is approximately .80 volts which is comparable with the 

32 

known value of Eg. 

Figures 8-12 and 8-13 give the C-V curves taken at 
two different gates on the sample. The major changes in 
these two curves from the first gate are: 1) the percentage 
change between accumulation and inversion capacitance is 
different (in Fig. 8-7 the change is 66%, in Fig. 8-12 
it is 50% and in Fig. 8-13 the change is 30%), 2) the 
minimum low-frequency capacitance value does not occur at 
the same voltage on all three curves, and 3) different 
values of accumulation capacitance. Possible explanations 
for these differences are nonuniformity in both the semi- 
conductor and insulator and, in the low-frequency case, 
errors in test measurements . The maximum value of capaci- 
tance is due to the different gate areas of the Hg probe. 

The differences among the three C-V curves is sig- 
nificant and, as a result, the data analysis curves differ. 
Figure 8-l4 gives the (C^/Cg)^ - Ug plot for the gate 3 
location. This figure and Fig. 8-8 look similar but are 
different because of the different scaling present. Two 
other differences exist between the data analysis curves 
of gate 3 and the two other gates. These differences are 
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For valid results, in equilibrium, the range of Vg 
should be equal to the value of the energy gap. For Ge 
Eg = 0.67V. 



130 




131 




Figure 8-13 • Experimental C-V Curves for Ge-MIS (p-Type) 
Gate #3 
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Figure 8-14. (C^^/Cg)^ - Ug Plot for Ge-MIS (p-Type) Gate #3 
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indicated in Figs. 8-15 and 8-l6. The flatband voltage has 
gone positive with a value of approximately 0.8V and 
begins to slope upward at a ip value of .28V. The reason for 
these differences is due to the shift of the voltage to the 
right for the C-V curves in Fig. 8-13. Nonuniformities in 
the insulator and semiconductor are probable causes of these 
differences . 

Table VII gives a listing of the tabular data output 
from the program. It can be seen from Table VII that there 



Quantity 




Gate #1 


Gate #2 


Gate 


HF Only 


ENS(X10^^m“ 


3) 


- 2.52 


- 2.98 


- 19.7 


— 


ENB(X10^^m“ 


3) 


-1.67 


- 3.56 


- 36.2 


— 


^b(v) 




-0.229 


-0.248 


- 0.308 


— 


N„„(X 10 ^^V" 
(at ip = 


V=) 

0) 


3.96 


13.6 


79.3 


2.73 


N^„„(X 10 ^^m 
"" (at Vg 


= 0 ) 


2.17 


1.71 


2.46 




Table VII. 


Summary of Tabular 
ent Gate Locations 


Output Data for Three 
on Ge-MIS (p-Type) 


Differ- 



are slight differences between gates 1 and 2 but there is a 
large difference in gate 3. This again suggests a nonunifor- 
mity in the semiconductor and insulator under gate 3. 

The interface state density data obtained from the 
high-frequency only C-V data analysis is also listed in the 
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table. The density obtained from this analysis agrees with 
that obtained above. 

2 . Analysis of n-Type Ge-MIS 

Based on the high-frequency only C-V data, results 
for the n-type Ge-MIS sample appear in Table VIII. Differ- 
ent gate areas led to similar differences as those described 
for the p-type Ge-MIS sample and will not be presented here. 
Table VIII shows that the interface state density is 5 
times higher than that for the p-type. The theoretical nor- 
malized high-frequency minimum capacitance is lower by an 
order of magnitude than the measured value. Figure 8-5 
shows that the experimental high-frequency curve had not 
reached a minimum. Higher bias voltage was not tried to 
avoid breakdown of the insulating layer. 
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Quantity 



Value 



0.303 ev 

L. 9.57X10"^cm 

a 





*^MINLF 


0 


C\J 

0 




^MINHF 


0 


.021 


Bias 


Surface Potential 


Interface State 


Energy Level 


(V.) 


(V.) 


Density (cm“2) 


of States (ev . ) 


-9.0 


0.24 


13 

2.50x10 


0.540 


-5.0 


0.25 


13 

1.43X10 


0.550 


-1.0 


-0.17 


12 

2.12X10 


0.130 


1.0 


-0.19 


12 

3.41X10 


0.110 


5.0 


-0.25 


1.48C10*^ 


0.050 


10.0 


-0.30 


13 

2.94X10 


0.003 


Table 


VIII . Summary of 
(n-Type) 


Tabular Output Data 


for Ge-MIS 
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IX. EFFECT OF ELECTRON BOMBARDMENT ON SILICON 



N-CHANNEL DEPLETION MOSFET 



A. INTRODUCTION 

The motivation for this section is twofold. The first 
factor is to test the computer program described in Chapter 
VI on a commercially tested standard Si-MOS. The second 
factor is to Investigate the effects of electron bombardment, 
which simulates the radiation environment around Jupiter 
[IM], on a Si-MOSFET. 



B. SAMPLES 

Commercial Si N-channel depletion MOSFET 's were selected. 
The 3.nd C-V^ characteristics were measured with a 

Tektronix type 576 curve tracer and Boonton capacitance 
meter respectively. The results of these measurements are 
shown in Figs. 9-1 and 9-2. 



C. EFFECT OF ELECTRON BOMBARDMENT 
1 . Electron Bombardment 

The bombardments were conducted at the Naval Post- 
graduate School Linear Electron Accelerator (LINAC). The 
energy of the electrons was 64.5 MEV. The electron beam was 



3 3 

The high-frequency C-V curve was obtained with the 
source and drain connected to ground, the substrate to the Hi 
terminal and the gate to the Lo terminal. The low-frequency 
curve was obtained with the source and substrate to the Hi 
terminal, the drain open, and the gate to the Lo terminal. 
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a) Positive Gate Voltage 
(bottom most curve: = O.V) 




b) Negative Gate Voltage 
(topmost curve: V^ = O.V) 

PlEure 9-1. Ijo-V.o Characteristics for MOSFET before 

^ aS dS Irradiation 
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Figure 9-2. C-V Characteristics for MOSPET before Irradiation 





I 




focused to provide a nearly circular pattern of approximately 
1 cm^ on a test target as viewed from remote TV cameras. 

The concentration or electron dosage was determined from the 
charge accumulated in a Faraday cup. This charge was stored 
in a .ly farad capacitor and the concentration found from the 
equation N = CV/q electrons per cm^, where V is the voltage 
across the capacitor. A secondary method to determine the 
total dose is to time the LINAC "on time" since, by experi- 
ence, a known dose exists in each pulse for a particular 
power setting. Electron dose rates from 10^^ e/cm^ to 
10^® e/cm^ in steps of a power of 10 increase were attempted. 
These dosages were not realizeable since difficulty was en- 
countered in the LINAC. Table IX gives the actual electron 
concentration obtained. 



Step 


Dosage per step (e/cm^) 


Total accumulated 






dosage (e/cm^ ) 


1 


1X10^^ 


1X10^^ 


2 


. 2X10 


3X10^^ 


3 


6oxio‘^ 


6.3X10*® 


4 


2.X10‘® 


2 . 063 x 10 *® 


Table IX. 


MOSFET e-Beam Total Accumulated 


Dosage Obtained 



2 . Experimental Results 

Two electrical characteristics of the MOSFETS were 
measured before and after the bombardment. These 
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curves. The 



characteristics are the C-AL and the I 

Lr 




pre -bombardment characteristics were given in section B. 

The MOSFET characteristics after bombardment steps 2, 3, and 
and 4, are shown in Figs. 9-3 through 9-8. Figures 9-3 



these figures that electron bombardment is, in effect, clos- 
ing the channel. For each electron dose step the drain-to- 
source current is lowered, for a certain gate voltage, and 
the valley of the low-frequency C-V curves widens. A second 
MOSFET was also tested with the first. The results for it 
are similar and will not be presented in this thesis. 

3 . Experimental Analysis 



of the MOSFET was obtained from the characteristics. 

It was found that the was shifted to the right after 



electron dosage is shown in Table X. 

Second, the physical changes of the device are pre- 
sented by determining the change of parameters of the MOS 
structure. These parameters are; 1) 'i’ 3 -VQ^ the surface 
potential and gate voltage relation, 2) the interface 



-jli 

Vth is the gate voltage at which the MOSFET begins to 
conduct . 



through 9-6 give the 




characteristics and Figs. 9-6 



through 9-8 give the C-V^ characteristics . It is clear from 



Two types of analysis were made. First the effect 

34 

of electron bombardment on the threshold voltage, 



electron bombardment. The shift of V, 



th 



as a function of 




a) Positive Gate Voltage 
(bottom most curve: V_ = O.V) 

(JT 




b) Negative Gate Voltage 
(topmost curve: V^ = O.V) 

MOSFET IdS“"'^dS Characteristics after a Total 
Accumulated Dosage of 3X10^^ e/cm^ 



mu 



Figure 9-3. 






a) Positive Gate Voltage 
(bottom most curve: = O.V) 




b) Negative Gate Voltage 
(topmost curve: V~ = O.V) 



Figure 9-4- MOSFET IdS“^dS Characteristics after a Total 
Accum.ulated Dosage of 6.3X10^^ e/cm^ 
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Positive Gate Voltage 



Figure 9-5* 



MOSFET IdS“^dS Characteristics after 
Accumulated Dosage of 2.063X10^® e/cm 
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Figure 9-6. MOSPET C-V Characteristics after a Total Accumulated Dosage 
of 3X10^^ e/cm^ 
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Figure 9 - 7 . MOSFET C-V Characteristics after a Total Accumulated Dosage 
of 6.3x10*^ e/cm^ 
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Figure 9 - 8 . MOSPET C-V Characteristics after a Total Accumulated Dosage 
of 2.063x10^® e/cm^ 



Dosage (e/cm^) 


hh 


0 


-l.il 5 


1X10 


-1.30 


2X10 


-1.20 


6.3X10 


-0.80 


2.063x10 


> 0 



Table X. MOSPET Threshold Voltage Shift as a Function of 
Total Accumulated Dosage 



state density and 3) effective charge in the insu- 

lator as a function of the surface potential. 

The relation between surface potential and applied 
gate voltage changed considerably. Figure 9-9 shows this 
relationship before irradiation and Pig. 9-10 shows the 
effect after the third radiation step. These two figures 
reveal that the flatband voltage becomes more negative a : i 
the entire curve shifts upward with irradiation. This cl inge 
is caused by the ionization traps created in the insulat r. 

The change in the interface state density versus ^ 
curve is represented in Fig. 9-11 before and 9-12 after the 
third radiation step. The trend after steps 1 and 2 was to 
shift the peak of the curve to the left. In Fig. 9-12, 
however, the peak which occurred at about -0.8V in Fig. 9-11 
has disappeared and another peak with an order of magnitude 
less at 2.5V has appeared. The disappearance of the first 
peak is unexplained. This peak was expected to continue to 
move to the left on increased irradiation. 
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Figure 9 - 9 * MOSFET V„-V Characteristics after a Total 
Accumulated^Dosage of 3X10^^ e/cm^ 
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The effective insulator charge also changes consid- 
erably after the third bombardment step. Initially 
was approximately 1X10 ^coul/m^ at zero surface potential 
and decreased gradually to 0.35V where it decreased sharply 

-3 

After bombardment was 2.5X10 at zero surface poten- 

tial decreased rapidly to O.IOV where it leveled off. 

4 . Discussion 

The effect of electron bombardment shows that the 
MOSPET is vulnerable to radiation damage to such an extent 
that it changed from a depletion to an enhancement MOSPET 
of poor quality. This drastic change occurred after an 
accumulated dosage of 2X10^^ e/cm^ which was received after 
an e-beam exposure time of 10 minutes. 

The electron bom.bardment closed the channel since 
more positive threshold voltages were required to support 
drain-to-source current as the dosage increased. This fact 
is apparent since the drain-to-source current for a parti- 
cular drain-to-source voltage at various gate voltages 
decreased as the dosage increased. 

The bombardment widened the valley in the low-fre- 
quency C-V curve and stretched out the accumulation to in- 

35 

version region on the high-frequency curve. The effect 
is caused by the creation of interface states and on ioniza 
tion traps in the insulator. With increased dosages the 



35 

Similar effects were produced by gamma radiation of 
0.5 to 0.7 megarads exposure to MOS capacitors [28]. 
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surface potential versus applied voltage shifted upward and 
the corresponding flatband voltage became more negative. 
Both the surface state density and effective insulator 
charge were affected considerably by the irradiation. 
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X. CONCLUSION 



Three MIS studies were made. 

(1) MIS of IV-VI compounds and alloy semiconductors: 

The purpose was to develop a fabrication procedure for 
MIS of PbTe, PbSe, PbSnTe and PbSnSe and to investigate 
whether their Interface qualities are adequate for infrared 
CCI applications. 

To fabricate MIS samples, 3 to lOy thick thln-films 
of these semiconductors were first deposited on CaP^ or BaF^ 
substrates followed by e-gun deposition of 100-450A® thick 
AljOj or SiOj insulating layers. Hg probes were used as 
metal gates. 

The C-V and C-t of the MIS samples were measured mostly 
at 300°K and on selected samples at 77°K. It was found 
that both AljOj and Si 02 layers showed considerable leakage. 
Capacitance measurements were distorted. In spite of this 
difficulty, C-V variations have been obtained which indicated 
that accumulation, depletion and inversion behaviors occurred 
and MIS of these IV-VI compounds and alloy semiconductors 
behave qualitatively like that of Si-MOS. 

However, quantitative comparisons of the measured C-V 
results with theoretical calculations revealed that the 
following dlscrepencles exist: 
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a. The measured change of MIS capacitance from accumu- 
lation to inversion is up to 5 times larger than calculated 
values in some cases. 

b. The measured voltage bias range of the transition 
from accumulation to inversion is up to 4 times larger than 
the calculated values. 

c. Some unexpected temperature variations of the MIS 
capacitance were found. It vjas found that the accumulation 
capacitance varied with the temperature. 

It should be noted that similar problems existed in the 
early stage of the development of Si-MOS. These discrepan- 
cies suggested that the following improvements of this 
research should be made; 

a. A deposition procedure of good insulator must first 
be developed. Insulators other that AI2O3 and SiOj having 
large dielectric constants should be used. It is believed 
that the leakage is the major cause of much of the 
difficulty . 

b. If good insulation materials can not be developed 
for use with these semiconductors, then the possible use of 
a double junction MIS should be investigated. 

c. After the successful development of good insulator 
or double junction MIS, the surface preparation of the semi- 
conductor and post treatment of MIS such as annealing can 
then be pursued to improve the quality of MIS for possible 
IRCCI applications. 
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d. Low-frequency measurement and lov;-temperature meas- 
urement must be Improved in order that the computer program 
developed in this thesis can be successfully used to analyze 
the basic physical parameters of these MIS which are neces- 
sary to build up a good understanding of their behaviors. 

(2) MIS of Ge: 

Both n and p-type Ge-MIS have been made by e-gun 
deposition of l80 to 250A° of AI 2 O 3 . Their resistivity are 
AOQ -cm and 0.05^ -cm respectively. C-V and C-t measure- 
ments Indicated that they behave qualitatively like Si-MOS. 
Analysis of these data showed that the interface state den- 
sities and the flatband voltages are: 

Quantity 0.05^ -cm (p) ^O.Ofl -cm (n) 

Nss(V"^m“^) 4X10^® 2. 5x10 

(V) -1.20 



The Ge-MIS C-V characteristics were not similar at all 
gate locations. It is suggested that sample fabrication 
techniques be improved to eliminate these nonuniformities. 

(3) Effect of electron bombardment on n-channel Si 
depletion MOSFET: 

65 .^Mev electron bombardments were carried 
out in the Naval Postgraduate School Linear Electron Acceler- 
ator. Doses of 10^^ j 6X10^^ , and 2X10*® were used. It was 
found that the most sensitive change is the decrease of 
threshold voltage. Before bombardment, the was -1.^5V 

after a dose of 2X10^* e/cm^ , V^^ became positive, which 
changed the MOSFET from depletion to enhancement. 
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U.U. 

<xc 



OCCCCOOCOCCCCOCO 



r,QAPH J\ Tl Fn 

GFRHANIIIM NA-I.1FI7 T»300 X«250A 

Knsr P-TVPE SUBSTRATE 

HAS BFFN PI rJTTED. 
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PROGRi^M TO COMPUTE THE SURFACE POTE NTI A L , I MPURI TY 
CONCENTRAT ION, INTERFACE STATE DENSITY, AND EFFECTIVE 
SURFACE CHARGE IM MIS STRUCTURES FROM READ IN 
EXPERIMENTAL CAPACITANCE VS VOLTAGE HIGH AND LOW 
FREQUENCY CURVES 



THE PARAMETERS USED INCLUDE(ALL MKS UNITS): 

AREA=AREA OF GATE 

BOLTZ = BOLTZMANN' S CONSTANT 

EPO=PERMITTI VITY OF FREE SPACE 

Q=ELECTRCNIC CHARGE 

ENI = I NTRINSIC CARRIER C ONS E NT RAT I ON 
T=TEMPERATURE 

ERS=DIELECTRI C CONSTANT OF SEMICONDUCTOR 
ERO=DIELFCTRIC CONSTANT OF INSULATOR 
XO=INSULATOR THICKNESS 

CCL=LOW FREQUENCY ACCUMULATION CAPACITANCE 
CLONEW=LOW FREQUENCY CORRECTED COL FOR AL203 
COH=HIGH FREQUENCY ACCUMULATION CAPACITANCE 
CLO=LOW FREQUENCY CAPACITANCE 
VA=LOW FREQUENCY GATE VOLTAGE 
CHI=HIGH FREQUENCY CAPACITANCE 
VH=HIGH FREQUENCY GATE VOLTAGE 

VS=POTENTIAL AT INTERFACE SURFACE OF SEMICONDUCTOR 
ENS=IMPURITY CONCENTRATION AT SEMICONDUCTOR SURFACE 
ENB=IMPURITY CONCENTRATION IN SEMICONDUCTOR BULK 
PF=PICOFARAD 

PHIB=POTENTI AL IN SEMICONDUCTOR BULK 
NSS=INTERFACE STATE DENSITY 
QEFF=EFFECTIVE INTERSTATE CHARGE 



DIMENSION VS (500 ) , VA( 500) , CH( 500 ) , VH( 500 ) , I D( 19) , 
1CL(500) ,CL0(500) , CHI (500) 

DATA IDN/4HIDN /,IDP/4HIDP / 

N=0 

1 READ(5,9,END=100) N AM E, ARE A , CLON EW 
9 FCRMAT(A4,t>X,E10.3,F10.3) 

IF(AREA.EQ.O. ) GO TO 1 
READ(5,10) B0LTZ,£P0,0,ENI 

10 F0RMAT(4E15.5) 

READ(5,11) T,ERS 

11 FORMAT( 2F10.3) 

VTH=( BCLTZ''T)/Q 
EPS=ERS4EP0 
READ(5,20) COL 

20 FORMAT ( FIO. 3) 

DO 30 I =1 , 500 
READ(5,25) VA(I),CLO(I) 

25 FORMAT (2F10. 3) 

IF(VA(I ). EQ. 500.0) GO TO 31 
CL(I)=CLD(I)/COL 

30 CONTINUE 

31 NL = I-l 
READ! 5, 40) COH 

40 FORMAT ( FIO .3 ) 

DO 50 1=1,500 
READ(5,45) VH( I) ,CHI( I) 

45 FORMAT(2F10.3 ) 

IF( VH( I ) .EQ.500. 0) GO TO 51 
CH( I ) =CHI ( I )/COH 

50 CONTINUE 

51 NH = I-l 

IF (AREA.GT. l.OE-10) GO TO 53 
READ(5,52) ERO,XO,PF 

52 FORMAT ( FIO. 3, 2515. 3) 

AREA=(C0HTPFfX3)/ ( ERCtEPO) 

WRITE(6,54) AREA 
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54 FORMAK 6X, E9.3) 

53 I F(NL . EQ.O .GR.NH. EQ.O .OR . I\1L .GT . 500.0R.MH.GT .500) GO 
ITO 1 
N = N+1 

I FCCLONEW.NE.O . ) COL = CLONEW 
READ(5,60) WHAT 
60 FORMAT ( FIO . 1 ) 

WHAT INDICATES TYPE OF PLOT DESIRED 
IF( WHAT.LT.O) POINT PLOT 

IF(WHAT .EQ.O ) POINT AND SMOOTH CURVE PLOT 
IF( WHAT.GT.O) SMOOTH CURVE PLOT 

CALL VSOFVAI VS, EN S , ENB , C OL , C L , VA , NL , C OH , CH , VH,NH, AREA 
1 , EPS, VT H, ENI, WHAT ) 

IF(NAME.EQ.ION.OR. NAME.EQ.I DP) GO TO 1 
CALL STATES! ENB, V S , CL , CH , VA , NL , ARE A , E P S , VTH ,ENI ,WHAT) 
CALL QEFF( ENS , EN B , VS , V A , NL , CH , C OL , AR EA, E P S , VTH, EN I 
1, WHAT) 

GO TO 1 
END 

100 CONTINUE 

SUBROUTINE VSOFVA PROVIDES PLOTS OF: 

(C/CS)v-"‘2 - US CLF - VL CHF - VL VS - VA 
AND PROVIDES TABULAR DATA ON: ENS ENB PHIB 



SUBROUTINE VSOF VA { VS , ENS , ENB , CO LO , CLN, VL , NL ,COHI , CHN , 
1VH,NH, AREA, EPS , VT H, EN I , WHAT ) 

DIME NS I ON VS (500) , CLN ( 5 00 ) , VL ( 5 00 ) , CHN ( 50 0 ) , V H( 500 ) , 
IF (500), G( 500),H( 500) ,X( 500) ,RANGE(4) 

DATA PI /3. 14159/, Q/ 1 .602E- 19/ , RUP/0. 75/ , RLO/0.95/, 
lPF/1. OE-12/ 

REALMS TITLK 12) ,TITL2( 12) , TI TL3( 12) 

REAL LABL1/'C-US'/,LABL2/'SFITV,LABL3/' LF '/ 

REAL LABL4/' HF ' / , LAB L5/ ' DE PL ' / , L ABL6/ ' R EGN* / 

REAL LABL7/'VSVA'/ 

READ(5,101) TITLl ,TITL2,TITL3 
101 FORMAT! 6A 8) 

CHNMIN=SMALL (CHN, NH) 

CUPPER = l.-RUPi'- (1 .-CHNMIN) 

CL0WER=1 .-RLO^I 1. -CHNMIN) 

NUS = 0 

CSMIN=CHNMIN<'COHI / (I. -CHNMIN) 

DO 1 1=1, NL 
F(I)=1 .-CLN( I) 

IF(F( I ).LT.O.) F(I ) = 0. 

1 CONTINUE 

WE NOW HAVE l.-C/CQ(LF) IN 
CALL SUMUP(F , VL, VS,NL) 

DO 2 1=1, NL 

F(I) =FTP( VL( I) , CHN,VH,NH) 

G(I) = ( l./F( I )-l. )-^2 

IF(F( I ) .GT.CUPPER .OR.F( I ) .LT.CLCWEP ) GO TO 2 

NUS=NUS+1 

H(NUS)=G( I) 

X(NUS)=VS(I ) 

2 CONTINUE 

F NOW CONTAINS C/CO(HF) WRT VL 

G NOW CONTAINS (C0/CS)»^*'2 WITH RESPECT TO VL OR VS 
H AND X CONTAIN POINTS FROM G AND VS FOR STLFIT 
VA = FTP(X(1),VL,VS,NL) 

VB = FTP( X(NUS) , VL,VS,NL ) 

CALL STLFIT( H, X,NUS, SLOPE ,YCEPT) 

USO=-YCEPT/SLOPE 
H( 2) =1 .lt-BIGST( 3 ,NL) 

X(2)=( h(2)-YCEPT )/SLOPE 
H(l)=-.2->H(2) 

X( 1) =( H( 1 )-YCEPT) /SLOPE 
IF(X( 1 ) .LT .X ( 2) ) GO TO 5 
XTR = X(l) 

X(l) = X(2) 

X(2) = XTR 
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XTR = H(l) 

HO.) = H(2) 

H(2) = XTR 

H AND X CONTAIN NOW THE STLFIT STRAIGHT LINE 
5 WRITEI6, 12) 

12 FORMAT Cl',' (C0/CS)»i-^2 - US PLOT •) 

RANGE! 1)=0.8 

RANGE(2 )=-0. 10 
RANGEO )=2 .0 
RANGE! 4)=-0. 1 
IF !WHAT) 50»51»52 

50 CALL UTPLOT!VS,G»NL, RANGE, 1,1) 

CALL UTPL0T(X,H,2,RANGE,1,3) 

GO TO 53 

51 CALL UTPLOT!VS,G,NL , RANGE, 1 ,1 ) 

CALL UTPL0T(X,H, 2,RANGE,1,3) 

52 CALL DRAW!NL,VS,G, 1,0,LABL1,TITL1,0.0,0.0,0,0,0,0, 

16, 9,1, LAST) 

CALL DRAW!2,X,H, 3,0,LABL2,TITL1,0.0,0.0,0,0,0,0,6,9, 
ILAST) 

53 CONTINUE 

DO 3 1 = 1, NL 

CLN! I ) = C0L0*CLN! I ) 

3 CHN! I ) = COHITFd ) 

CLN AND CHN NOW CONTAIN C!LF) AND C!HF) WRT VL 
WRITE!6,13) 

13 FORMAT! CLN G CHN - VL PLOT ') 

F!l) = FTPIVA, CHN,VL,NL ) 

F!2) = . 75FCHNMIN*C0HI 

GO) = VA 

G!2) = VA+l.E-3 

H!l) = FTP!VB,CHN,VL,NL) 

H!2) = F!2) 

X!l) = VB 
X!2) = VB + l.E-3 

F,G,H,ANDX NOW CONTAIN DEPL. REGION MARKERS 

RANGE ! 1 ) =4.0 

RANGE! 2 )=-8.0 

RANGEI3 ) = C0L0 + 1.0 

RANGE (4 ) = .65LCHNMIN4-COHI 

IF !WHAT) 60,61,62 

60 CALL UTPLCT(VL,CLN,NL, RANGE, 1, 1) 

CALL UTPL0T(VL,CHN,NL,RANGE,1,2 ) 

CALL UTPLOT(G,F,2,RANGE,l,2) 

CALL UTPL0T!X,H,2,RANGE,1,3) 

GO TO 63 

61 CALL UTPLOTIVL, CLN, NL, RANGE, 1,1) 

CALL UTPLOT!VL, CHN, NL , RANGE, 1, 2) 

CALL UTPLOTIG ,r ,2 ,RANGE ,1 ,2 ) 

CALL UTPLOT!X,H, 2,RANGE, 1, 3) 

62 CALL DRAW(NL,VL,CLN,1,0,LABL3,T ITL2,4.0, 0.0,0,0,0,0, 

16, 9,1, LAST) 

CALL DR AW!NL, VL, CHN,2, 0,LABL4,TITL2,4.0, 0. 0,0,0,0,0, 

16. 9.1, LAST) 

CALL DRAW!2,G,F, 2,0,LABL5,TITL2 ,4.0,0. 0,0,0,0,0, 

16. 9.1, LAST ) 

CALL DRAW! 2,X,H ,3,0,LABL6,TITL2 ,4. 0,0. 0,0, 0,0,0, 

16.9.1, LAST) 

63 CONTINUE 

ENS=-2.’"COHO^‘2/! EPS^QiARE Af' i«2T SLOPE) T P F^'^2 
ENB = SIGN!QEY!CSMIN-^PF, AREA, VTH, EPS, ENI ) ,ENS) 

PHIB = SIGN! VTUStALOG! ABS ! ENB)/ ENI ), ENS ) 

WRITE! 6, 89) 

89 FORMAT! ' 1' , ' VS - VA • ) 

WRITE16,90) ENS, ENB, PHIB 

90 FORMAT!//,' I MP UR I T Y C ONC ENTR A T I ON ' , / / 3 X , ' E NS = ', 
11PE9.2,' PER CUBIC MET ER ' , / 8X , • ENB = ',E9.2, 

2' PER CUBIC METER* ,/8X,' PHI B = ',0PF6.3,' VOLTS') 

VS0 = -VTH1>ABS! ENB ) / ENS- ! ENS-ENB ) /ENS^PHI B/2. 

DELTA=VSC-USO 

DO 4 I=1,NL 

4 VS!I)=VS! I )+DELTA 
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RANGE (1)=4.0 
RANGE( 2)=-8.0 
RANGE( 3 )=1 .0 
RANGE (4 l=-0. 5 
IF (WHAT) 70,71,72 

70 CALL JTPLOKVL, VS, NL, RANGE, 1,0) 

GO TO 73 

71 CALL UTPLOT(VL, VS, NL, RANGE , 1,0) 

72 CALL DRAW(NL, VL, VS ,0,0 ,LABL7,T ITL3, 0.0, 0. 0, C, 0, 0, 0, 
16, 9,1, LAST) 

73 CONTINUE 
100 RETURN 

END 

SUBROUTINE STATES PROVIDES A PLOT OF NSS-PSI AND 
TABULAR DATA ON NSS AT THREE VALUES OF PS I 



SUBROUTINE STATES ( ENB, VS , CL F, CHF, VA, NLF , AREA, EP S , VTH, 
lENI ,WHAT) 

DIMENSION VS (500 ) ,CLF( 500) , CHF( 500) , VA( 500) ,PSI ( 500) , 
1ENSS(500) , RANGE(4 ) 

DATA PI /3. 14159/ , PF/1. OE-1 2 / , Q/ 1 . 602 E- 1 9 / 

REALMS TITL4(12) 

REAL LABL8/' NSS’ / 

READ( 5, 104) TI TL4 
104 FORMAT (6A8) 

COL=BIGST(CLF,NLF ) 

COH=BIGST(CHF,NLF) 

NSS=0 

F=PF/( AREA'-O) 

PHIB = VTH^^AL0G( ABS(ENB)/EN1 i 
IF(ENB .LT.O. ) PHI B=-PHI6 
U1=AMIN1(-1. 7-PHI 3,0. ) 

U2=AMAX1(-1.7"PHIB,0. ) 

DO 10 1=1, NLF 

IF(VS( I ) .LT.U1.0R.VS( I ).GT.U2)G0 TO 10 
NSS=NSS+1 

PSI(NSS)=-(PHIB+VS( I) ) 

ENSS(NSS ) = (CDLi-CLF( I )/(C0L-CLF ( I ) ) -CCH> CHF (I ) / 

1 (COH-CHF( I ) ) )#F 
IF(ENSS(NSS).LT.O.)ENSS(NSS)=0. 

10 CONTINUE 
IF(NSS.GT.2) GO TO 5 
WRITE! 6,90) 

90 FORMAT!///’ VS DATA INADEQUATE FOR NSS CALCULATION') 
GO TO 100 
5 LAST = NSS/2 
DO 11 1=1, LAST 
J = NSS-I+1 
X = PSI ( J) 

PSI ( J ) = PSI ( I) 

PSK I ) = X 
X = ENSS! J) 

ENSS ( J ) = ENSS( I ) 

11 ENSS( I ) = X 

X = FTP( 0. ,ENSS ,PSI ,NSS) 

Y=0. 

WRITE (6 ,8) 

8 FORMAT! ’1’,' NSS - PSI PLOT') 

W RITE ! 6,9 ) ENSS ! 1), PSI ! 1),X, Y, ENSS! NSS), PSI (NSS) 

9 FORMAT!//,’ INTERFACE STATE DENS ITY ’ / / 5 X , 3 HNS S , 10 X , 

13HPSI/6X, ’! /VMv‘J'2 ) ! VOLTS ) ’ //3 ! 6 X , 1 P£9 . 2 , 

22X,0PF8 .3/ ) ) 

RANGE! 1 )=0.3 
RANGE!2 ) = -0.2 
RANGE!3 )=4 .0E17 
RANGE! 4)=4.0E15 
IF (WHAT) 80,81,82 

80 CALL UTPLOT!PSI , ENSS, NSS, RANGE, 1 ,0) 

GO TO 83 

81 CALL UTPLOT!PSI, ENSS, NSS, RANGE, 1,0) 
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82 CALL DRAW(NSSiPSI ,ENSS,0,0,LABL8,TITL4,0.1,0.0,0,0,0,0 
1,6,9, 1,LAST) 

83 CONTINUE 
100 RETURN 

END 

C SUBROUTINE OEFF PROVIDES A PLOT OF QEFF-VS AND TABULAR 

C DATA ON QEFF NEFF AT THREE VALUES OF VS 



SUBROUTINE QEFF( ENS , ENB, VS , VA , N LF , CHF , COL , ARE A, E PS , 
1VTH,ENI ,WHAT) 

DIMENSION VS (530 ), VA(500) ,CHF( 500),QFF( 500) ,VF( 500) , 
1RAN6E( 4) 

DATA Q/1.602E-19/ ,PI/3. 14159/, PF/l.OE-12/ 

INTEGER QEFF/4H0EFF/, L0G0/4HL0GQ/ 

REAL«^8 TITL5( 12) 

REAL LABL9/'QEFF '/ 

REA0(5,105) TITL5 
105 FORMAT! 6A8) 

COH=BIGST (CHF, NLF ) 

F = Q§:EPS»^AREA«‘'^2 
R=2.=t ( ENI/6NB )1s1^2 
NFF=0 

PHIB = SIGN(VTHf'ALOG(ABS(ENB)/ENI ) ,ENB) 

U1 = AMIN1(0.,-1.7'^PHIB) 

U2=AMAX1(0.,-1 .7-PPHIB) 

DO 10 1=1, NLF 

IF(VS( I ) .LT .Ul.OR.VS! I ) .GT.U2) GO T3 10 

NFF=NFF+1 

VF(NFF ) = VS( I ) 

X=VF(NFF)/VTH 
IF(ENB.LT.C. ) X=-X 

QFF(NFF) = (COL^( VS( I )-VA( I ) ) v PF-F* ( COH-CHF ( I ) ) / (Coni' 
1CHF( I)t'PF)"M ENS-ENB'!'(EXP(X )+RtS INH(X) ) ) )/AREA 
10 CONTINUE 

IF(NFF.GT.2) G3 TO 5 
WRITE(6 ,90 ) 

90 FORMAT! ///• VS DATA INADEQUATE FOR QEFF CALCULATION') 
GO TO 100 
5 Y=-PHIB 

X=FTP( Y,QFF ,VF ,NFF) 

S1 = ABS (QFF( 1 )/Q ) 

S2=ABS (X/Q) 

S3=ABS (OFF(NFF) /Q ) 

WRITE(6,8 ) 

8 FORMAT! '1',' QEF!= - VS PLOT') 

WRITE(5,9)QFF( 1 ) , S 1 , VF ( 1 ) , X , S2 . Y , QFF ( F ) , S3 , VF { NF F ) 

9 FORMAT!///' EFFECTIVE SURFACE 'CHARGE, STATE DENSITY* 
1//8X, 4HQEFF , 8X, 4HNEFF , 6X, 2HVS/5 X, • (COUL/M^^'2) ’ ,3X , 

2' !/M4*2 ) • , 3X, ' (VOLTS) '//3!6X, 1PE9.2,1X,E11.2, 
31X,0PF8.3/) ) 

RANGE! I )=0.5 
RANGE (2 )=-0. 1 
RANGE ! 3 )=1.0E-2 
RANGE!4) = -1.0E-3 
IF (WHAT) 85,85,87 

85 CALL UTPLQT!VF,QFF ,NFF, RANGE, 1,0) 

GO TO 88 

86 CALL UTPLOT(VF, QFF,NFF, RANGE, 1, 0 ) 

87 CALL DRAW!NFF,VF,QFF,0,0,LABL9,TITL5,0.0,0,0,0,0,0,0, 
16,9, 1, LAST ) 

88 CONTINUE 
100 RETURN 

END 

C SMALL FINDS THE SMALLEST OF X 
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FUNCTION SMALL(X,N) 

DIMENSION X(500) 

SMALL = X (1 ) 

DO 1 I=2,N 

IF(X( I ) .LT.SMALU SMALL = X(I) 

1 CONTINUE 
RETURN 
END 

C SUMUP INTEGRATES BY TRAPEZOIDAL RULE 



SUBROUTINE SUMUP ( Y, X , F, NPT ) 

DIMENSION Y(500) ,X(500) ,F(500) 

F(l) = 0. 

DC 10 1=2, NPT 

lOF(I) = FU-1) + ( X( I )-X( I-l ) )*(Y ( I )+Y( I-in/2. 
RETURN 
END 

C FTP SETS THE MAXIMUM UPPER AND LOWER LIMITS ON Y 



FUNCTION FTP(VAL,Y,X,NPT ) 

DIMENSION X( 500) ,Y(500) 

DATA IC/1/ 

IF(IC.GT.NPT) IC=1 
IF( VAL-X( IC) ) 1,2,3 

2 FTP = Y(IC) 

GO TO 5 

3 IF( IC.EQ.NPT) GO TO 2 
IC = IC+1 

IF(VAL-X( I C) ) 6,2,3 
1 IF(IC.EQ.l) GO TO 2 
6 IC = IC-1 

IF( VAL-X( IC) ) 1,2,4 

4 FTP = Y(IC) + ( VAL-XdOIS'IYdC + D-YdC) ) /(X(IC+1)-X(I 

5 RETURN 
END 

C STLFIT FITS A STRAIGHT. LINE TO A CURVE 



SUBROUTINE STLFIT( Y, X,N, SLOPE , YCEPT) 
DIMENSION Y(500 ) , X(500) 

XBAR = EXPECT(X,N) 

Q = 0. 

A = 0. 

DO 20 1=1, N 
DEV = Xd )-XBAR 
Q = Q+DEV*^2 
20 A = A + DEVAY( I ) 

IFIQ.NE.O.) SLOPE = A/Q 
IF(Q.EQ.O.) SLOPE = 0. 

YCEPT = EXPECT(Y,N) - SLOPE^^XBAR 

RETURN 

END 

C EXPECT PROVIDES THE SUMMATION OF Z 



FUNCTION EXPECT(Z,N) 
DIMENSION Z(500) 

EXPECT = 0. 

DO 10 1=1, N 

10 EXPECT = EXPECT+Z(I) 
EXPECT = EXPECT/FLOAT (N) 
RETURN 
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END 



C BIGST FINOS THE LARGEST OF X 



FUNCTION BIGST(X,N) 

DIMENSION X(500) 

BIGST = X(l) 

DO 1 1=2, N 

IF(X( I ) .GT. BIGST) BIGST = X(i) 

1 CONTINUE 

RETURN 

END 

C QEY SOLVES FOR BULK IMPURITY CONCENTRATION 



FUNCTION QEYCSMI N , A REA , VTH , EPS , ENI ) 

DATA P 1/3. 14159/, Q/1.602E- 19/.PF/1.0E-12/ 
VAL=VTH*(CSMIN/AR£A)i<*2/( EPSAQ^ENl) 

N = 0 
X = 13. 

DX=1. 

1 X=X+DX 

IF (EXP (X/2. )/SQRT (2.-t(X-2. ) )-VAL ) 1, 100, 2 

2 DX=+.li'DX 
N=N + 1 

IFIN.GT .2 ) GO TO 100 

3 X=X-DX 

IF(X.LT.3.) GO TO 100 

I F( EXP (X/2. )/SQRT (2.^(X-2. ) )-VAL) 4, 100, 3 

4 DX=.1=J^DX 
GO TO 1 

RETURN 

100 QEY=2.«-VTHfc(X-2. )*(CSMlN/AREA)A*t2/( EPS=?Q) 
END 
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PROGR^^M TO CAL:uLATE THE MINIMUM HIGH FREQUENCY 
CAPACITANCE OF A 4-6 SEMICONDUCTOR, WITH 
VALUES OF SEMICONDUCTOR DIELECTRIC CONSTANT, 
VERSUS IMPURITY CONCENTRATION 



THE INPUTS TO THIS PROGRAM ARE 

resc=semicondu:tor relative dielectric constant 

REINS=INSULATOR RELATIVE DIELECTRIC CONSTANT 
T=INSULAT0R thickness in ANGSTROMS 
DA=IMPURITY DOPING LEVEL 



REALt^^S TITLE 

REAL^^S LABEL(4),TITEL (10, 12),TITLEI( 12) ,LABLE 
DIMENSION DA(20) , DL0GDA(20) ,CMINH(20) 

N=13 

READ(5,2) (DA( I), 1=1, N) 

2 FORMAT! 10E6. 1) 

READ! 5, 30) (LABEL! I ), 1 = 1,4) 

30 FORMAT (4A8) 

DO 100 1=1,5 

READ!5, 31 ) !T ITEL! I ,M) ,M=1, 12) 

31 F0RMAT!6A8) 

DO 20 L=l,12 

T ITLEI ! L)=TITEL! I,L) 

20 CONTINUE 

READ!5,1) REINS, T, TITLE 
1 FORMAT (2F10. 2, A8 ) 

WRITE (6,21 ) 

21 FORMAT! •!' ) 

WRITE!6,7) TITLE, T 

7 FORMAT!' I NSULA TOR ' , A8 , • WI TH A THICKNESS 0F',F7.0, 

1 • ANGSTROMS ' ) 

WRITE(6,6) REINS 

6 FORMAT! 'OINSULATOR RELATIVE DIELECTRIC CONST ANT F7 . 1 ) 



DO 95 J=l,4 
REA0!5,3) RESC 

3 FORMAT ! FIO .2 ) 

WRITE!6,8) RESC 

8 FORMAT!//,' SEMICONDUCTOR RELATIVE DIELECTRIC CONSTAN* 
1 ,'T' ,F7.1,/) 

WRITE! 6,4) 

4 FORMAT (4X, 'DOPING' ,/, 5X, 'LEVEL ’ ,2X, 'C! MIN )HF' ,11X, 

1 ' XMAX' ,8X , ' EG' ,3X, ' 2 *^DELEV' , // ) 

IFU-3) 11,11,12 

11 MC=2 

IF(J.EO.l) MC=1 
GO TO 13 

12 MC=3 

13 DO 90 K=1,N 
LABLE=LABEL! J ) 

TEMP=300. 

VTH=8. 625E-05feTEMP 
ES = RESC=«-8.85E-14 
X=0.0 

EG=0.181+4.52E-04*TEMP-0. 5 68F X+ 5 . 8^ X*:?4 
EMD=1 . 12'i-EG 

ENV=4. 82E15^ ! (TEMP^ EMD)**! .5 ) 

EINS=RE INS^-8. 85E-14 
DELEV = VTH^'ALQG( ENV/DA!K ) ) 

DELEV2 = 2.1= DELEV 

XMAX= SORT! 2. i^E S/1. 6E- 19/DA !K)J^( EG-2. «'DE LE V) ) 

CSMIN=ES/XMAX 

CINS=1.0E08^-EINS/T 

CMINHIK )= l./( l.+CINS/CSMIN) 

DLOGDA! K)=ALOG10(0A(K) )-16.0 
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WRITE ( 6,5) DA(K) ,CMINH(K) ,XMAX, EG,DELEV2 
5 FORMAT (1PE10.1,0PF10.4,E15..4, 2F10.4) 

90 CONTINUE 

CALL DRAW( N,DLG3DA,CM1 NH, MC, J, LA BLE, TITLE!, 0.5, 0.0,0, 
10, 1,2, 9, 6,0, LAST) 

95 CONTINUE 
100 CONTINUE 
END 
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